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^ ■ Abstract 

CN \ We consider the possibility that the hghtest supersymmetric particle is a 

■ heavy gluino. After discussing models in which this is the case, we demon- 

\^ . strate that the 5-LSP could evade cosmological and other constraints by 

\ virtue of having a very small relic density. We then consider how neutral and 

' charged hadrons containing a gluino will behave in a detector, demonstrating 

00 ' that there is generally substantial apparent missing momentum associated 

\ with a produced ^-LSP. We next investigate limits on the ^-LSP deriving 

' from LEP, LEP2 and RunI Tevatron experimental searches for excess events 

in the jets plus missing momentum channel and for stable heavily-ionizing 
Qh. charged particles. The range of m~ that can be excluded depends upon 

(— I \ the path length of the g in the detector, the amount of energy it deposits 

in each hadronic collision, and the probability for the g to fragment to a 
pseudo-stable charged hadron after a given hadronic collision. We explore 
\ how the range of excluded m~ depends upon these ingredients, concluding 

' that for non-extreme cases the range 3 GeV < m~ < 130 — 150 GeV can 

be excluded at 95% CL based on currently available OPAL and CDF anal- 
yses. We find that Runll at the Tevatron can extend the excluded region 
(or discover the g) up to m~ ~ 160 — 180 GeV. For completeness, we also 
analyze the case where the g is the NLSP (as possible in gauge-mediated 
supersymmetry breaking) decaying via g ^ g + gravitino. We find that the 
Tevatron RunI data excludes m~ < 240 GeV. Finally, we discuss application 
of the procedures developed for the heavy ^-LSP to searches for other stable 
strongly interacting particles, such as a stable heavy quark. 

1 Introduction 

In the conventional minimal supergravity (mSUGRA) and minimal gauge-mediated 
(mGMSB) supersymmetry models, the gaugino masses Mj at low energy are pro- 
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portional to the corresponding and are in the ratio 

M3 : M2 : Ml ~ as : 02 : ai , (1) 

as would, for example, apply if the Mj evolve to a common value mi/2 at the GUT 
scale Mjj in the SUGRA model context. However, well-motivated models exist in 
which the Mj do not obey Eq. (0). In particular, the focus of this paper will be on 
models in which M3 is the smallest of the gaugino masses, implying that the gluino 
will be the lightest supersymmetric particle (LSP). (We note that we explicitly do 
not consider g masses as low as those appropriate in the light gluino scenario [0], 
which some [0] would claim has now been ruled out.) 

One such model is the O-II string model in the limit where supersymmetry 
breaking is dominated by the universal 'size' modulus P, ^ (as opposed to the 
dilaton). Indeed, the O-II model is unique among the models considered in in 
that it is the only string model in which the limit of zero dilaton supersymmetry 
breaking is consistent with the absence of charge/color breaking. In the absence 
of dilaton supersymmetry breaking, the gaugino masses arise at one-loop and are 
therefore determined by the standard renormalization group equation coefficients 
and by the Green Schwartz parameter 5gs- The O-II model in this limit results in 
the ratios 

Ms : M2 : Ml -(3 + 5gs) : (1 - Sgs) : (y - Sgs) , (2) 

and a heavy gluino is the LSP when 6gs ~ ^3 (a preferred range for the model). 

In the GMSB context, the possibility of a heavy ^-LSP has been stressed in 
Ref. 1^. There, the g is the LSP as a result of mixing between the Higgs fields and 
the messenger fields, both of which belong to 5 and 5 representations of SU(5), 
which are, in turn, contained in lO's of SO(IO). The basic idea is as follows. First, 
one implements the standard mechanism for splitting the color-triplet members of 
the Higgs from their SU(2)-doublet partners in the 5, 5 representations using an 
'auxiliary' 10. As a result of this splitting, the Higgs color triplets mix with the color 
triplet members of the auxiliary 10, both acquiring mass of order the unification 
scale, M[7. If one now identifies the fields in the auxiliary 10 with the messenger 
sector 10 fields, it is the messenger sector fields that supply the standard doublet- 
triplet Higgs splitting and whose color triplet members acquire mass ~ Mjj. As a 
result, the color-triplet messenger fields naturally become much heavier than their 
SU(2)-doublet counterparts. Since the masses of the gauginos arise in GMSB via 
loop graphs containing the messenger fields of appropriate quantum numbers, the 
result is that the (colored) gluino mass is suppressed by {M/MuY compared to the 
other gaugino masses, where M is the typical mass of a doublet messenger field. 
One requires that M/Mu ^ 0.1 in order to adequately suppress baryon number 
violating interactions mediated by the Higgs triplets (which are controlled by an 
effective mass of order M^/M). 
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Early outlines of the phenomenological constraints and possibilities for a heavy 
^-LSP appear in H, 0, H, ^ 0. Here, we attempt to refine these phenomenolog- 
ical discussions. For our phenomenological studies, we will make the assumption 
that all supersymmetric particles are substantially heavier than the ^-LSP.Q This 
is a conservative assumption in that discovery of supersymmetry will be easier in 
scenarios in which some of the other superparticles are not much heavier than the 
gluino. 

The outline of the rest of the paper is as follows. In section 2, we demon- 
strate the sensitivity of the relic gluino density to assumptions regarding the non- 
perturbative physics associated with gluino and gluino-bound-state annihilation. 
In section 3, we examine how energetic massive gluinos produced at an accelerator 
will be manifested in a typical detector. In section 4, we consider the constraints 
from LEP and LEP2 data on a massive gluino produced in e~^e~ — > qqgg. In section 
5, we examine constraints on a massive ^-LSP from the existing RunI data in the 
jets plus missing momentum channel and explore the prospects for improvements 
at Runll. In both sections 4 and 5, we discuss how the constraints/limits depend 
on the manner in which a ^ is manifested in a detector. We consider limits on a 
heavy ^-LSP that arise from searches for heavy stable charged particles at OPAL 
and CDF in sections 6 and 7, respectively. In section 8, we present Tevatron limits 
on a gluino that is the NLSP of a gauge-mediated supersymmetry breaking model, 
decaying via g g + gravitino. In section 9, we outline possible applications of 
the procedures developed for the heavy gluino to other new particle searches, in 
particular searches for a stable heavy quark. Section 10 presents our conclusions. 
The reader is encouraged to begin by scanning the concluding section 10 so as to 
get an overview of our results and the issues upon which he/she should focus while 
working through each section. 

2 The Relic Gluino Density 

Before embarking on our discussion of direct accelerator limits, it is important to 
determine if a massive gluino LSP can have a relic density that is sufficiently small 
to be consistent with all constraints. In particular, as discussed in [U, 13; H S 0' 
its relic density must be sufficiently small that it cannot constitute a significant 
fraction of the dark matter halo density. Otherwise, it would almost certainly have 
been seen in anomalous matter searches, underground detector experiments and so 
forth. We will show that non-perturbative physics can lead to large enhancements 
in the relevant annihilation cross sections, with the result that the relic density 
could be very small. 

We begin with a very brief review of the standard approach for computing a 
relic density. First, one determines the freeze-out temperature Tp, which is roughly 

^This is natural for the sfermfons in the O-II model, since the uiq SUSY-breaking scalar mass 
parameter is automatically much larger than mi/2- 
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the temperature at which the annihilation rate for two gluinos falls below the rate 
at which the universe is expanding. The standard form of the freeze-out condition 

Here, Gn is Newton's constant, x = m~/T, (7^ = 2 x 8 is the number of gluino 
degrees of freedom, and g*{T) is the density degree-of-freedom counting factor. In 
all our computations, we employ the exact formula of Ref. |11] for (cr^'^^f): 



1 r°° 

9 y g 



where /? = y^l — 4m~/ s is the velocity of the g^s in the initial state center-of-mass 

frame; (a^^'^v) is computed numerically. The above [a'^'^^v) form assumes only that 
the ^'s (or i?°'s, see below) remain in kinetic equilibrium for all temperatures (as 
seems highly likely given that they re-scatter strongly on either quarks/gluons or 
hadrons, respectively, even after freeze-out). We then numerically integrate the 
Boltzmann equation. Defining as usual Y = n~/s (where s is the entropy density 
and n~ is the gluino number density), the standard result is 
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(5) 



where the subscript (F) refers to current (freeze-out) temperature and h*(T) is 
the entropy degree-of-freedom counting factor .0 As usual, l/Yp <^ l/^o and can 
be neglected. Finally, we compute the current gluino mass density as 



po = m-^o = m~soYo = m~/i*(To)— Tq Yq , (6) 
and 

Pc 8.0992 X 10~47 GeV^ ' ^ ' 

The estimates in the literature H, ^, |[| for the relic density of a massive gluino 
differ very substantially, at least in part due to different assumptions regarding the 
size of the annihilation cross section. Perturbatively, the annihilation cross section 
is o-f,""" = a{gg gg) + Eg (^{M ^ ql) with: 

<^{gg - gg) = ^ { log ^ 21 - 60' - 3(3' - 33(3 + 17(3' 




a[QQ ^ QQ] = —^[6- p-)[6- IJ ). (9) 

lops 



^Note that only standard model particles are counted in computing g* and h* since all super- 
symmetric particle are presumed to be heavier than the 5. 
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[In Eq. ® , /3 = y 1 — 4mg/s, rriq being the quark mass.] We observe that as /5 — ^ 0, 
/3(T™" approaches a constant unless the employed is allowed to increase in a 
non-perturbative manner. (Note that this is in sharp contrast to the /3(T™ oc 
p-wave behavior for the XiXi annihilation cross sections; since the ggg vertex does 
not contain a 75, gg annihilation can occur in an s-wave and is much stronger at 
low p.) For our perturbative computations we employ (Q) evaluated at Q = y/s, 
where a^(Q) is the usual moving coupling, oc 1/ log((5^/A^) at one loop. (When 
employed at small Q, see below, (y^{Q) = 1 will be the maximum value allowed.) 

However, near the threshold, ^/s ~ 2m~, non-perturbative effects can be ex- 
pected to enter. There are many possibilities. Consider first multiple gluon ex- 
changes between interacting ^'s. These will give rise to a Sommerfeld enhancement 



factor |12, 13, 14, 15|, which we will denote hj E, as well as logarithmic enhance- 



ments due to soft radiation |jl5[. Here, we retain only E, which takes the formQ 

-1 



P 



1 — exp 



Cnag 



(10) 



with C being a pro cess- dependent constant. The Egg (Egg) for gg gg {gg qq) 
is given by taking C = 1/2 (C = 3/2). If one examines the derivation of E, 
then one finds that the typical momentum transfer of the soft gluon exchanges 
responsible for £" is Q ~ P^g- Thus, we choose to evaluate E using as(/3m~).0 
The C values quoted above are those appropriate to color averaging in the initial gg 
state. Color averaging is relevant since the high scattering rate of gluinos (off gluons 
etc.), continually changes the color state of any given gluino, and, in particular, 
does not allow for the long time scales needed for the Sommerfeld enhancement to 
distort 113 the momenta of the relic gluinos so that they become organized into 



color-singlet pairs with low relative velocity. In what follows, we will employ the 
shorthand notation Eap^'^ = Egga{gg — >• gg) + Eqga{gg — >• qq). 

As an aside, we note that multiple soft-gluon interactions between the final 
state q and q in gg ^ qq result in a repulsive Sommerfeld factor at small f3 (since 
the qq are in a color octet state). However, this is not an important effect since the 
gg —>■ qq cross section vanishes as /5 — > anyway. We do not include this final-state 
Sommerfeld factor in our calculations. 

We will consider two possibilities for computing Eap^'^. In the first case, ap'^^ 
is computed using 0^(1/5) and E is computed using a^{Pm~), with the result that 
/3£'(T™" oc 1//? as /5 — i> 0, recalling that a^{j3m~) is not allowed to exceed 1. In 
the second case, we employ a 'non-perturbative' form for a^, denoted a^^, defined 

^The Sommerfeld enhancement factor takes the form l + CTras/{2f3) for small Cnas/{2(3). We 
extend this to the region of large CTTas/{2(3) by using the standard exponentiated form given. 
In the perturbative next-to- leading order results of [|5|, Er^ and Eqg are evaluated at the 
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factorization scale fi. In the perturbative expansion approach, a next-to-next-to-leading order 
calculation is required to determine the appropriate effective scale at which to evaluate the next- 
to-leading Sommerfeld factor. 
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by replacing 1/ log((5^/A^) in the form by l/log(l + Q^/A^). (This form is 
that which corresponds to a roughly linear potential a large distance, and was first 
discussed in Ref. [|TB| with regard to the charmonium bound state spectrum.) cTp""^ 
and E are evaluated using a^^{y/s) and af^(/5m~), respectively. The result is that 
(5Ea^^ oc at small (3. In both cases, the growth of Ea^'^ will be cutoff by 
requiring that Ea^'^ not exceed Ea^"^ = l3~^/m^, the largest annihilation cross 
section that we wish to consider. 

Of course, as is well-known from the charmonium analogue [T^, the Sommerfeld 
enhancement at best provides an average (in the dual sense) over the resonance 
structure that is likely to be present. Further, just as in charmonium, the Som- 
merfeld enhancement is a precursor to the formation of gg bound states that will 
occur once the temperature falls below the typical binding energy. This binding 
energy would be of order ~ ct'irn- to the extent that short range Coulomb-like 
color attraction is most important, but terms in the potential between the two 
gluinos (that possibly rise linearly with the separation) can also play an important 
role. Thus, it is difficult to be precise about the temperature at which this tran- 
sition occurs, but it is almost certainly above the temperature of the quark-gluon 
deconfinement transition. If gg bound state formation were to be complete, the an- 
nihilation rate n-a^^'^v (where n~ is the number density of gluinos per unit volume) 
would be replaced by the decay rate for the gg bound state. In the charmonium 
analogy, this decay rate is proportional to |A^p|\&(0)p, where Ai is the matrix 
element for the decay, oc a^/m|, and |^I/(0)| is the magnitude of the wave function 

at the origin, oc [agm-]^/^. The result is a decay rate proportional to a^m-. The 
important feature of this result is that the bound state draws the two gluinos to- 
gether (as represented in |\I^(0)|) so as to overcome the perturbative behavior of the 
annihilation |A^p. A full treatment would have to implement a coupled-channel 
treatment in which the gg bound state formation would be treated in analogy with 
the standard approach to e~p recombination in the early universe. Those ^'s that 
are not absorbed by gg bound state formation prior to the temperature falling 
below the deconfinement transition temperature would end up inside bound states 
containing one g and one or more gluons or light quarks; most likely the = gg 
bound state would be dominant. The rate of annihilation of the i?°'s is far from 
certain (as discussed below). Although we [1^ are exploring the possibility of 
implementing this full scenario, there are clearly many uncertain ingredients. We 
presume that the resulting relic density will be bracketed on the high side by the 
Sommerfeld enhancement result and on the low side by the limit where very few 
gg bound states form before the confinement transition, below which strong R^R'^ 
annihilation takes over. 

In this latter extreme non-perturbative scenario, we imagine that at small f3 
there will be a transition where the ^'s condense into color singlet bound states 
containing one g and light quarks and/or gluon(s); as noted above, we shall assume 
here that the lightest is the = gg. (An electrically charged LSP bound state has 
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much stronger cosmo logical constraints and is easier to see at accelerators.) For j3 
above the transition point, we will employ a^'^^ without any enhancement factor 
E. For /3 below the transition, the appropriate annihilation cross section will be 
that for R^R^ n's. It is often assumed (see, e.g., H, |^, H) that the non- 
perturbative cr™ will be = A(3~^ /m^, where the [3~^ factor is the standard 
result for s-wave annihilation of spin-0 particles and A is an uncertain constant not 
too different from unity. We will consider this possibility even though we regard 
such a large annihilation cross section as being unlikely since annihilation must 
remove the gluino quanta, implying, in a parton picture, gluino exchange in the 
t-channel.0 Note that if A scales as l/m-, we would obtain ~ cTp""^ (both 
behaving as as /? — and having similar normalization); the result would 
be a relatively smooth transition as the temperature crosses the deconfinement 
boundary, yielding a result not very different from our perturbative case (with no 
Sommerfeld enhancement factor). 

In our numerical work, the choice of = A(3~^ /m^^ with A = 1 is labelled as 
I. As an alternative, we also consider a second choice (II): a'^^ = I/itl"^, such that 
/?CT™ vanishes (like Up"'^) as /3 — >■ 0. Although II has no particular model moti- 
vation (other than representing a kind of average of s-wave and p-wave behavior) , 
it allows us to assess the importance of the small /? behavior of crl^. We will see 
that it leads to significantly larger relic densities than I. For a given choice of cr^p^, 
the exact point of transition between ap^^ and and its smoothness are also 
crucial ingredients in determining the relic density. 

• For the transition point we consider two choices: (a) the total gg kinetic 
energy in the center-of-mass falling below a given limit L, with L ~ 0.2 — 
1 GeV (we employ L = 1 GeV in our numerical results — the relic density 
increases with decreasing L); (b) twice the g momentum falling below L. We 
note that the transition occurs roughly at /5 ~ ^ L/m~ and {3 ~ L/m~ in 
cases (a) and (b), respectively. To the extent that the condensation of ^'s 
into bound states is controlled by the typical temperature, the KE criterion is 
the most natural. It is because it leads to large increases in the rehc density 
that we have considered the more moderate (b) possibility. 

• For the smoothness of the transition we also consider two options: (i) use ap'^^ 
for larger (3 with an abrupt transition to the non-perturbative annihilation 
form for f3 below the appropriate limit; (ii) a smooth transition in which 
(Tp""^ is evaluated using a^^{Q) and Q is taken to be the net kinetic energy, 
i/i — 2m~, or 2p^^ in cases (a) and (b) above, respectively. The modified a^"^ 
is employed until it exceeds cr^, after which point the latter is employed. A 

^In, for example, the model of Ref. for strong scattering, A would scale as for 

annihilation, in sharp contrast to the R'^R^ /J^i?" scattering cross section which would scale 
with the inverse size squared of the gg bound state (which would have comparable size to a pion 
or proton bound state). 
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smooth transition will lead to larger relic density than the sudden transition 
choice. 

Altogether, we shall consider eleven cases. The first three are: (1) cr™" {E = 1); 
(2) Edf;"'' with E as given in Eq. evaluated using a^{Q = /3m~); and (3) Eap'''' 
with E computed using (y^^{Q = f3m~)] in (2) and (3) E'er™" is not allowed to 
exceed Ea^"^ = (3~^ /m^. The remaining eight cases are specified by various 
scenarios: (4) (I,a,i); (5) (II,a,i); (6) (I,b,i); (7) (II,b,i); (8) (I,a,ii); (9) (II,a,ii); (10) 
(I,b,ii); (11) (II,b,ii). 



0.040 




m~ (GeV) 



Figure 1: Tp/m~ as a function of m~ for the 11 cases described in the text. The 
solid lines correspond to results for cases (1), (2) and (3), respectively, in order 
of decreasing Tp. Results for cases (4) (I,a,i) (5) (II, a, i) (6) (I,b,i) (7) (II,b,i) are 
the lower dashed, dotted, dot-dot-dash and dash-dot lines, respectively. Results 
for cases (8) (I,a,ii) (9) (II,a,ii) (10) (I,b,ii) and (11) (II,b,ii) are the upper dashed, 
dotted, dot-dot-dash and dash-dot fines, respectively. Tfiis figure assumes L = 
1 GeV, see text. 

Results for the freeze-out temperature and the relic gluino density for the eleven 
cases detailed above are shown in Figs. |I] and |^, respectively. As expected, the 
freeze-out temperature for a relic gluino (relative to the mass m~ of the gluino 
relic) is lower (by roughly a factor of two) than in the case of a weakly-interacting 
relic particle. The ordering of the curves for the eleven different cases can be easily 
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m~ (GeV) 



Figure 2: Qh"^ as a function of m~ for the 11 cases described in the text. Line 
notation as in Fig. |I], with sohd hues for cases (1), (2) and (3) in order of decreasing 

understood on the basis of the strength of the annihilation cross section for each 
case as a function of (3. 

After freeze-out takes place, annihilation remains substantial (especially in cases 
where a^'^^ jumps to a large value at small f3) and the relic-density continues to 
decline. The current relic density is thus very strongly dependent upon the model 
employed. Fig. |^ shows that Qh"^ can be substantial (even corresponding to an over- 
closed universe for m~ > 10 TeV) if a purely perturbative approach is followed, or it 
can be extremely small out to very large m~, as in case (I,a,i) where = P~^/ml 
and an abrupt transition from a^'^ to based on the KE criterion is employed.^ 
Almost any result in between is also possible. Further, the second sub-electroweak 
scale inflation discussed by some (see, for example, Ref. |2^) would dilute even 
the purely perturbative relic densities to an unobservable level. Until the non- 
perturbative physics issues can be clarified, and late time second inflation can be 
ruled out, we must assume that the relic g (or more properly R^) density is small 
enough that constraints from anomalous nuclei in seawater, signals associated with 

^This and the other related cr™ cases evade the upper bound on the mass of the dark matter 
particle of Ref. , based on s-wave dominance of the cross section and partial wave unitarity, 
by virtue of the fact that cr^p' ^ /3~^/mJ ^ Trf3~'^/m^ (the latter being the s-wave unitarity 

limit) can arise from, for example, the coherent contribution of many partial waves. 
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annihilation in the core of the sun, interactions in underground detectors etc. are 
not significant. In the following sections, we discuss the extent to which accelerator 
experiments can place definitive constraints on the heavy ^-LSP scenario. 

3 How a heavy gluino LSP is manifested in de- 
tectors 

Before turning to accelerator constraints on the ^-LSP scenario, we must determine 
how a stable gluino will manifest itself inside a detector. This is a rather comph- 

cated subject. The important question is how much momentum will be assigned 
to the jet created by the gluino as it traverses a given detector. This depends 
on many ingredients, including, in particular, the probability P that the gluino 
fragments to a charged i?-hadron, R^, vs. a neutral i?-hadron, It is useful to 
keep in mind the following two extremes. 

• Very little energy would be assigned to the g if it always fragments into an 

which interacts only a few times in the detector and deposits little energy 
at each interaction. 

• Large energy would be assigned to the g if it undergoes many hadronic inter- 
actions as it passes through the detector, with large energy deposit at each 
interaction, and/or if it fragments often to a R"^ following a hadronic collision. 
In particular, when the g moves with low velocity through the detector while 
contained within an R^, it will deposit a substantial fraction of its energy 
in the form of ionization as it passes through the calorimeters. Further, for 
non-compensating calorimeters this ionization energy is overestimated when 
the calorimeter is calibrated to give correct energies for electrons and pions. 

In addition, in the OPAL analysis to be considered later, if the gluino R- 
hadron is charged in the tracker and at appropriate further out points in 
the detector, it will pass cuts that cause it to be identified as a muon, in 
which case the momentum as measured in the tracker is added to the energy 
measurement from the calorimeter and a (much too small) minimal ionization 
energy deposit is subtracted from the calorimeter response. In this case, the 
energy assigned to the g 'jet' can actually exceed its true momentum. 

In all our discussions, it should be kept in mind that in current analysis procedures 
jets or jets containing a muon are always assumed to have a small mass, so that 
the momentum of a jet is presumed to be nearly equal to its measured energy. 

3.1 Hadronic energy losses: the g ^ R'^ case 

In this subsection, we explore the energy loss experienced by a heavy g passing 
through a detector as a result of hadronic collisions. An early discussion of the 
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Figure 3: Average energy loss, (AE), in a collision as a function of f3 for the three 
cases described in the text. Results are shown for m/jo = 5, 25 and 140 GeV. At 
high f3, curves are ordered according to increasing m^o. 



issues appears in Ref. [PT|. These would be the only energy losses if the g almost 



always moves through the detector as part of an state. (This would be the 
case if charge-exchange reactions are significantly suppressed because the charged 
g bound states are substantially heavier than the or if the states undergo 
rapid decay to an state.) The first question is how much energy will the R^ 
lose in each hadronic collision as a function of its current (3 value. As a function 
of |t| and (where t is the usual momentum transfer invariant for the R^ and 
mx is the mass of the system produced in the R^N — > R^X collision) the energy 
loss is given by 

^j,^rnl-m% + \t\ 

2mN ^ ' 

where we have assumed that the appropriate target is a single nucleon rather 
than the nucleus as a whole or a parton (both of which are estimated to be irrelevant 
in l2ll]). To estimate the average /\E per collision, we must assume a form for 
We have examined three different possibilities: 



d\t\dmx 



;i) ^^1^ oc 1 for \t\ < 1 GeV^ and zero for |t| > 1 GeV^. 



(2) ^i^i'^'^j given by a triple-Pomeron form |22 



X 



ap(0)-l 



where ap(|t|) = 1 - 0.3|t| and (3{\t\) = 1/(1 + |t|/0.5 GeV^)^ is a typical 



parameterization. For the parameterization of Eq. ([TJ), the result for the av- 
erage energy loss {AE) is independent of the maximum value (if > 0.5 GeV^) 
allowed for \t\. 
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(3) ^^g^ oc 1 for |t| < 4 GeV" and zero for \t\ > 4 GeV. 

We compute the average value of AE as a function of the (3 of the in the rest 
frame of the target nucleon: 

where |t|mm,max('^x) = 2[E{mN)E{mx) Tp{mN)p{mx) - m\o\ with E{m) = {s + 
mj^o — m?)/{2y/s) and p{m) = A^/^(s, m^o, m^)/(2v^) [with X{a,b,c) = + 6^ + 
— 2(a6 + ac + 6c)], where s = m|,o + + 27mKom7v [with 7 = (1 — In 
integrating down to mx = ttin in Eq. (0), we include both elastic and inelastic 
scattering (using the same cross section form).[] We note that the above kinematic 
limits for |t| as a function of mx must be carefully incorporated in order to get 
correct results for (AE); in particular, |t|mm I^Uax as mx a/s — mRO. 



The results for {AE) obtained from Eq. (|13D in the above three cross section 
cases are plotted in Fig. ^ for three masses that will later prove to be of interest: 
mpto = 5, 25 and 140 GeV. We note that (AE) as a function of f3 is almost 
independent of the mass so long as m^o > 5 GeV. In what follows we will use 
the m^o = 25 GeV results for (AE) for all m^o. 

In order to understand whether any of the three models for ^^^J^^^ is reason- 
able, and, if so, which is the most reasonable, we examined the results given by our 
procedure in the case where the R^ is replaced by a pion. In so doing, the pion is 
viewed as retaining its identity (aside from possible charge exchange) as it traverses 
the detector, slowing down after each hadronic collision by an amount determined 
by the (AE) for the then current j3 of the pion. In our approach, since the elas- 
tic cross section is effectively included in our cross section parameterizations, the 
average distance between hadronic interactions of the pion is characterized by its 
path length At (in the notation of Ref. |23|) in iron (Fe) as determined by the total 
cross section. (We will also need to refer to the inelastic collision length, denoted 
by A/.) In Fig. |^, we show how the energy of a 100 GeV pion deteriorates to below 
5% of its initial energy as it undergoes successive hadronic collisions separated by 
At, using cross section models (1) and (2).0 In Fig. 24.2 of Ref. results for the 



number of A/ = 17 cm interaction lengths in iron required for 95% of the kinetic 
energy of a pion to be deposited as a result of hadronic collisions are given as a 
function of initial energy. We have computed this number for the (AE) predictions 
of our three cross section models; note that in our approach, hadronic interactions 
occur every At = 11 cm. The results^ for cross section models (1) and (2) are 

''For large /3 > 0.95, the purely elastic scattering component gives smaller (AE) than the 
inelastic scattering component. This should be incorporated in a more complete treatment. 

*Note that the (AE) values in Fig. ^ are not correct for a light hadron; we employ Eq. (|l^ ) 
computed numerically for the current (3 value just prior to a given collision. 

^Results are independent of whether the pion is assumed to be charged or not; i.e. dE/dx 
losses are not important. 
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Figure 4: We plot the energy of an incident 100 GeV pion after a certain number 
of hadronic collisions for the case (1) and (2) cross section models. 



given in Fig. |^ along with the results from Fig. 24.2 of Ref. For moderate 

energies, Fig. |^ shows that the triple-Pomeron case (2) yields rough agreement, 
but at higher energies predicts that 95% containment requires more A/ than ex- 
perimentally measured. The case (1) cross section predicts 95% containment for 
fewer A/ than actually measured for all initial energies. (Case (3) would predict 
that even fewer Xj would be required for 95% containment.) 

As we shall see, the main issue for detecting a ^-LSP signal is the amount 
of kinetic energy of the ^'s i?-hadron that is not deposited in the calorimeter. 
Deposited energy has many critical impacts in the context of the experimental 
analyses that we will later employ. We mention two here. First, for an event that 
is accepted by other cuts, larger missing kinetic energy implies a stronger missing 
momentum signal. This is the dominant effect for a ^-jet that propagates primarily 
as part of a neutral i?-hadron bound state. For the OPAL and CDF jets + missing 
momentum signals, considered in later sections, case (1) would then be conservative 
in that it leads to smaller missing momentum. Second, for larger missing kinetic 
energy a ^-jet that is propagating as a charged i?-hadron will be more frequently 
identified as being a muon. In the CDF jets + missing momentum analysis, muonic 
jets are discarded. As a result, case (2) will weaken this CDF signal for a charged 
i?-hadron (but not the jets + missing momentum OPAL signal, for which muonic 
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Figure 5: We plot the number of A/ = 17 cm {i.e. in iron) path lengths required 
for 95% containment of the energy of a pion. Experimental results from the PDG, 
Fig. 24.2 of Ref. are compared to predictions based on Eq. (|13[) for the case 
(1) and (2) cross section models. 



jets are retained). In later sections, we will use case (1) as part of our normal 
scenario-1, or "SCI", choices. Clearly, it will be important to explore sensitivity 
to the {AE) case choice. Of course, the net amount of energy deposited by a ^-jet 
is also influenced by the path length. At, of the g. As discussed below, a simple 
model suggests that A^ for the g is longer than A^ for a pion. For the graphs of this 
section, we will use the value of At = 19 cm derived from this model (see below). 
In later sections, however, we will discuss sensitivity to doubling and halving At 
relative to this "SCI" value. 

Turning now to the we compute the number of collisions Ncow required 
to deplete a certain percentage of the -R°'s initial kinetic energy. We carry out 
this computation by starting the out with a given (3 and stepwise reducing 
its kinetic energy according to the (AE) given in Fig. ^ Results for (AE) cases 
(1) and (2) are plotted in Fig. ^ for m^o = 5, 25 and 140 GeV. It is clear from 
this figure that what is important is how the initial (3 correlates with m^jo in the 
experimental situations of interest. The initial /3's that will be of relevance for 
these masses (which will prove to be of particular interest) are: (3 ~ 0.95 — 0.99 for 
fTT-RO ~ 5 GeV at LEP and m/jo ~ 25 GeV at the Tevatron; and f3 ~ 0.5 — 0.8 for 
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Figure 6: Number of collisions, A'^cou, required for an R'^ of the indicated mass (in 
GeV units) to deposit 90% or 50% of its kinetic energy given the initial (3 plotted 
on the X axis. The upper and lower hues of a given type are for {A.E) cases (2) 
and (1), respectively. The last P point plotted is P — 0.99. 



iriRO ~ 25 GeV at LEP and mf>o ~ 140 GeV at the Tevatron. In all cases, we see 
that a substantial number of collisions are required in order that the BP deposit a 
large fraction of its kinetic energy as a result of hadronic collisions. 

To interpret the above results it is necessary to know the number of hadronic 
coUisions that the is hkely to experience as it passes through the detector. Fur- 
ther, it is important to know how much of the energy deposited in a given hadronic 
collision will be measured as visible energy and, therefore, used in determining the 
energy of the associated 'jet'. In assessing the latter, we employ the following 
approximations. 

• For a neutral i?" (which interacts strongly only — no ionization), we presume 
that the energy deposited in both clastic and inelastic hadronic collisions in 
the calorimeters will contribute to 'visible' energy in much the same way as 
do energy losses by a pion. In this case, the calorimeter (which is calibrated 
using pion beams) will correctly register the amount of energy deposited 
by the R^. This should probably be more thoroughly studied in the case 
of elastic collisions for which all the energy deposited resides in recoiling 
nucleons which could have a somewhat different probability for escaping the 
absorbing material and creating visible energy in the scintillating material. 

• We assume that the energy deposited in uninstrumented iron, such as that 
which separates the calorimeters from the muon detection system in the CDF 
and DO detectors, is not visible. 
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For our cross section models, the number of hadronic collisions of the as it 
passes through the detector is determined by the total (and not just the inelastic) 
cross section for scattering on the detector material. This is normally rephrased 
in terms of the interaction length At in iron (Fe). The average number of collisions 
is then given by the number of equivalent Fe At interaction lengths that charac- 
terizes the detector. (However, it is conventional for detectors to be characterized 
in terms of their thickness expressed in terms of the number of inelastic collision 
lengths. A/, in Fe.) For the pion (which we take to be representative of a typical 
light hadron), we have already noted that Xt{'^) ~ H cm and A/(7r) ~ 17 cm 



The equivalent CDF and DO detector 'thicknesses' are specified in terms of the 
number of A/(7r). For all but a small angular region, the DO detector thickness 
ranges from 13 — 19 A7(7r), depending upon the angle (or rapidity) (the smallest 
number applying at rj = and the larger number at ~ 1.5). However, of this, a 
large fraction is in the CF or EF toroid magnets and is uninstrumented. The in- 
strumented thickness in which energy deposits are recorded ranges from ~ 7A/(7r) 
at ?7 = to ~ 9A/(7r) at r/ ~ 1.5). The CDF detector thickness at r] = consists 
of about 4.7A/(7r) of instrumented calorimetry and ~ 2.9A/(7r) of uninstrumented 
steel in front of the outer muon chamber. The instrumented portion of the muon 
detection system is fairly thin and will lead to little energy deposit. The LEP de- 
tectors have similar thickness for the instrumented category. In particular, at rj = 
OPAL has about 2A/(7r) of electromagnetic calorimetry and about 4.7A/(7r) in the 
instrumented iron return-yoke hadron calorimeter. Further, no additional unin- 
strumented iron is placed between the magnet return yoke and the muon detectors 
(which are drift chambers). To summarize, instrumented thicknesses at rj = are 
~ 5A7(7r) for CDF, ~ 6.5A/(7r) for OPAL and ~ 7A7(7r) for DO. At rj = 1.5 the 
thickness is perhaps as large as 9A7(7r) at DO. For r] < 1, uninstrumented sections 
add about 3A7(7r) for CDF and 6A/(7r) for DO in front of the muon chambers. To 
get the number of AT(7r) that corresponds to a given number of A/(7r), multiply 
the latter by ~ 1.6. Thus, the 5 (CDF), 6.5 (OPAL) and 7 (DO) A/(7r) for small 77 
convert to roughly 8 (CDF), 10 (OPAL) and 11 (DO) AT(vr). At ~ 1.5 add about 
3 AT(7r) to the CDF and DO numbers and perhaps 2 AT(vr) to the OPAL result. 
Uninstrumented thicknesses for 77 < 1 are ~ 5AT(vr) (CDF) and ~ 10AT(7r) (DO). 
OPAL has no additional uninstrumented iron prior to its muon chamber. 

We must now correct these thicknesses for the relative size of ajio ^ as compared 
to cr^AT, using the fact that Xt{tt) oc 1/crJ^. To estimate cr^o^v! "we employ the two- 
gluon exchange model for the total cross section developed in detail in Ref. [|TB|. 
Compared to the vrA^ cross section, the R'^N cross section must be increased by the 
ratio of Ca/Cf = 9/4 to account for the color octet nature of the i?" constituents, 
and it must be multiplied by (rj^o) / {r^), where (r^) is the (transverse) size-squared 
of the particle. In the simplest approach, which has substantial phenomenological 
support, (r^) is inversely proportional to the square of the reduced constituent 
mass of the bound state constituents: (r^) oc 4/mg vs. (rj^o) oc l/m^ (for m~ ^ 
mg), where and nig are constituent light quark and gluon masses, respectively. 
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Taking them to be similar in size, we find (T^otv ~ (9/16)crJ^, yielding At(-R°) ~ 
(16/9)AT(7r) ~ 19 cm. Using the factor of 9/16, and rounding up, the 8 (CDF), 
10 (OPAL) and 11 (DO) Xt{t^) instrumented thicknesses at small rj convert to 
5 (CDF), 6 (OPAL) and 7 (DO) Xt{R^). About 2 Xt{R^) should be added for 
rj ~ 1.5. For 1] < 1, about 3 (CDF) or 6 (DO) Xt{R^) uninstrumented interactions 
occur before the reaches the outer muon detection chambers. Below, we present 
results for 6, 7 and 8 instrumented hadronic interactions, as appropriate for the 
average measured energy deposit of i?°'s in the rj < 1.5 region at CDF, OPAL and 
DO, respectively. For later reference, it is important to note that the 8 hadronic 
interaction results are also appropriate for the total energy lost (even though not 
all is measured) due to hadronic collisions before reaching the outer (central) muon 
chambers at CDF. 

Obviously, a refined analysis by the detector collaborations to improve on the 
above will be quite worthwhile. More important, however, is understanding the 
extent to which the m~ region that can be excluded experimentally is sensitive to 
Xt{R^)- This will be examined when we consider exclusion limits based on OPAL 
and CDF analyses. 




Figure 7: The fraction of kinetic energy retained by the R^ is plotted as a function 
of its initial /3 for the cases of A^cou < 6, 7 and 8 for m^o = 5, 25 and 140 GeV. 
Upper and lower curves for a given mass are for (AE) cases (2) and (1), respectively. 

Our results for the fraction of the kinetic energy that is not deposited in 
the calorimeter (which will be the same as one minus the fraction included in the 
visible ^-jet energy /momentum) after A^^coii = 6, 7 and 8 hadronic collisions are 
presented in Fig. ^ as a function of the initial (3 of the Below, we make several 
observations that will be useful for understanding borderline cases that will arise 
in subsequent sections. 
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For OPAL at LEP (recalling that the number of hadronic collisions of the 
in the OPAL detector is close to 7): 

• For a 5 GeV with large f3 ~ 0.98, the triple-Pomeron [case (2)] (AE) 
implies that 7 interactions will deposit only about 20% of the kinetic 
energy. The constant cross section case (1) (AE) implies that about 45% of 
the KE would be deposited in 7 interactions. 

• For mRO = 25 GeV, A^cou = 7 and initial /3 > 0.5, the case (2) [(1)] cross 
section form would predict that no more than 20% [60%], respectively, of the 
i?° kinetic energy would be deposited in the calorimeter. 

For our CDF Tevatron analysis: 

• For m^o = 25 GeV and initial f3 > 0.95, less than 8% of the KE would be 
deposited in 6 interactions for the case (2) triple-Pomeron parameterization 
and less than 15% for the case (1) constant cross section choice. 

• For m^o = 140 GeV and initial f3 > 0.5, no more than 5% [10%] of the i?°'s 
KE would be deposited in case (2) [(1)] and contribute to visible energy in 
the detector. 

The key overall observation is that, in all large fraction of the gluino's 

kinetic energy will not contribute to visible energy in the detector. 

We now specify how events containing a stable R^ must be treated at the 
parton level in the standard OPAL and CDF analyses of the jets plus missing 
momentum channel that will be of special interest in what follows. The procedure 
given below assumes that the calorimeter calibration is such that energy deposited 
in the calorimeter by hadronic interactions is correctly measured. (This should 
be the case given that calorimeter calibration is established using a pion beam of 
known energy.) 

• As usual, in each event the visible three-momentum for a q, q or g jet is 
taken equal to its full three-momentum and its energy is taken equal to the 
magnitude of its three-momentum. 

• The visible energy of a ^ (as measured by the calorimeter) is taken equal 
to the total energy deposited in the instrumented calorimeter due to the ^'s 
hadronic collisions. 

• The magnitude of the three-momentum assigned to a 5? is taken equal to its 
visible energy {i.e. as if the visible ^-jet were massless) and the direction of 
the three-momentum is given by the direction of the g. 

• The invisible or missing momentum three-vector is computed as minus the 
vector sum of all the final-state three-momenta as defined above. Only trans- 
verse missing momentum is relevant for the experimental analyses. 
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• As usual, the absolute magnitude of the missing transverse momentum is 
termed the invisible or missing transverse energy. 

An alternative way of thinking about this is that for each ^-jet one computes the 
missing momentum as the difference 

IPtruel - IPapparentl = rn~{P-f - X{-f - 1)) , (14) 

where X is the fraction of the ^'s kinetic energy deposited and measured in the 
calorimeters of the detector: |papparent| = -'^ x KE — Xm~{'y — 1). The direction 
of a given ^'s contribution to missing momentum is the direction of the g. Note 
that even if X = 1, i.e. all the kinetic energy is seen by the detector, we find miss- 
ing momentum associated with the ^-jet of magnitude m~ (l — ^(1 — /?)/(! + , 
which is substantial for large m~ unless (3 is small. 

In the LEP and Tevatron analyses it will be important to note that since gf's 
are produced in pairs and in association with other jets with significant transverse 
momentum, the net missing momentum from combining the missing momenta of 
the two ^'s will not generally point in the direction of either of the ^-jets. Thus, 
^-pair events will normally pass cuts requiring an azimuthal or other separation 
between the direction of the missing momentum in the event and the directions of 
the various jets. 



3.2 Ionization energy deposits and the g — ^ R possibility 

We must now consider the possibihty that the g does not fragment just to an i?° 
that propagates through the detector without charge exchange. It might also have a 
significant probability for fragmenting to a (pseudo-stable) charged state, R^, when 
initially produced and after each subsequent hadronic interaction in the detector. 
(An example of an state would be a gud bound state.) We will assume that the 
initial and subsequent fragmentation probabilities are all the same. (We denote 
the common probabihty by P.) This would be the case if each time the i?-hadron 
containing the g undergoes a hadronic interaction in the detector the light quarks 
and/or gluon(s) are stripped away and the g then fragments independently of the 
previous i?-hadron state. A simple model for estimating P is the following. First, 
assume that the g is more likely to pick up a quark-antiquark pair to form a mesonic 
i?-hadron than three quarks to form a baryonic i?-hadron. If u, d {u, d, s) quark and 
antiquark types are equally probable, then of the 4 (9) possible quark-antiquark 
pairs only 2 (3) are charged and P = 1/2 (1/3) if the probability for fragmentation 
to gg is zero. Of course, if the = gg bound state is the lightest i?-hadron 
or is at least very close in mass to the gqq i?-hadrons, we expect that this latter 
probability is actually quite significant. If we assign the g a probability equivalent 
to all the quark-antiquark pair combinations included above, then P = 1/4 (1/6) 
in the u,d {u,d, s) cases, respectively. Thus, it would seem that P < 1/2 is quite 
likely. In considering the R^ states and the various neutral i?-hadron states on 
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a similar footing, we are implicitly assuming that all are stable against decay as 
they traverse the detector, i.e. that their lifetime is longer than ~ 10""^ sec. This 
will not be the case unless all the mass differences between the various states are 
smaller than mvr. Current estimates for the mass differences are too uncertain to 
reliably ascertain whether or not this is the case . 




Figure 8: The fraction of kinetic energy retained [i.e. that is not deposited) by a 
singly- charged g bound state is plotted as a function of its initial (3 for the cases 
of A^cou < 6, 7 and 8 for m = 5, 25 and 140 GeV. Upper and lower curves for a 
given mass are for (AE) cases (2) and (1), respectively. 

It is useful to consider first the extreme where P = 1 and compute the total 
amount of energy deposited, including both hadronic interactions and ionization. 
The hadronic energy losses are presumed to be the same as already discussed for the 
R^. For the ionization energy losses we employ the standard result for dE/dx from 
Ref. |2^. As before, we will parameterize the detector in terms of its equivalent size 



as if entirely made of Fe. Our procedure will be to integrate the ionization energy 
loss up to the point of the first hadronic collision at distance A^. The hadronic 
energy loss at this first collision will be computed for the then current (3 following 
our earlier procedures. We then integrate dE/dx starting from the f3 value retained 
by the after this first collision over a second of distance, compute the energy 
loss for this 2nd hadronic collision using the new current /3, and so forth. We will 
consider, as before, a certain number of hadronic collisions, A^cou = 6, 7 or 8. The 
At employed will be 19 cm, as discussed above. Ionization energy loss will be 
computed for Ncow segments of length At. The results corresponding to our earlier 
Fig. 1^ are presented in Fig. ^. There we plot, as a function of initial (3, and for 
-^cou = 6, 7 and 8, the fraction of kinetic energy of a singly-charged gluino bound 
state that is not deposited, after allowing for energy losses both from hadronic 
collisions and from ionization. 
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Figure 9: The minimum velocity /5min required for a singly- charged g bound state to 
retain non-zero kinetic energy after A^^coii = 7 or 8. The former (latter) is a rough 
estimate of what is required to penetrate to the OPAL (CDF) muon chambers. 
Results are plotted for (AE) cases (1) and (2). 




From Fig. |^ we see that for low enough /5 the will be stopped in the detec- 
tor. (For smaller initial jS, the ionization energy losses are larger and the velocity 
decreases rapidly.) This will be important when considering limits on a ^-LSP 
coming from searches for a stable charged particle that is heavily-ionizing. For ex- 
ample, CDF has placed strong constraints on such a stable charged object if its P 
is small enough for the particle to be at least twice minimal-ionizing (as measured 
soon after leaving the interaction vertex) but large enough that it will penetrate to 



the outer muon chamber p5[. For a singly-charged state, twice minimal-ionizing 
requires /37 < 0.85 or /?< 0.75. At CDF, roughly Ncoii = 8 collisions are ex- 
perienced by the charged hadron containing the gluino before reaching the outer 
central muon detector system. Fig. ^ shows that for m~ ~ 140 GeV (m~ ~ 25 GeV) 
P > 0.4 (> 0.6), respectively, is required in order that the g not lose all its kinetic 
energy before reaching the outer muon chamber. A plot as a function of m~ of 
the minimum initial /3, f3min, needed in order that the g retain non-zero KE after 
7 (8) collisions, and, therefore, penetrate to the OPAL (CDF) outer muon cham- 
bers, respectively, is presented in Fig. ^ Results are given for both the energy loss 
case (1) and case (2) models. We will later employ the lower limits for Ai'cou = 8 
and case (1) in assessing our ability to observe a charged gluino bound state as a 
penetrating heavily- ionizing particle in the Tevatron CDF experiment. 

Of course, if the g fragments part of the time to a neutral hadronic state 
and part of the time to a charged state and/or if charge exchange occurs as a 
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result of hadronic interactions, i.e. if P < 1 in the model discussed earlier, the 
results for energy loss and Prnin will be intermediate between the neutral and purely 
charged cases discussed above. However, in obtaining the accelerator limits based 
on heavily-ionizing tracks, to be discussed later, the reduced value of /3min that 
would apply for P < 1 is not important since the typical jS for the produced 
gluinos is substantially above P^in for the cases of interest. 



3.3 The momentum experimentally assigned to the ^^-jet: 
general g ^ R^, case 

Let us now return to the visible energy associated with P > probability for g 
appearance as an R^. In the case of a ^ traversing the detector and sometimes (or 
always) appearing as an R^, the procedure for determining this visible energy is 
analysis- and detector-dependent. 

First, we must note that both the OPAL and CDF hadronic calorimeters are 
constructed out of iron layers. These are intrinsically non-compensating in that 
purely ionization energy losses contribute more to the output energy measured by 
the calorimeter than do hadronic collision losses. For example, the CDF calorimeter 
is calibrated so that a 50 GeV pion beam is measured to have energy of 50 GeV. 
Using this same calibration, a 50 GeV muon beam is measured |^ to deposit 



2 GeV of energy whereas its actual energy loss as computed using the standard 
dE/dx of a muon in iron is only ~ 1.3 GeV. We define the ratio of calorimeter 
response to actual dE/dx loss from ionization as r. From the above, r = 1.6 for 
iron. The ionization energy deposited by an as it moves through the iron will 
be converted into r times as much measured calorimeter energy (which will be 
included in the visible energy /momentum of the ^-jet). The net energy deposited 
in the calorimeter after one complete interaction length will be measured to be 
-E^caiorimeter = ^Eionization + -Z^hadronic, after including the hadrouic energy deposit at 
the end. 

The next important consideration is whether there is a track, associated with 
the ^-jet, that is identified as a muon. 

• In the CDF jets + missing energy analysis discussed later, the ^-jet would be 
declared to be "muonic" if:|^ a) the g emerges from the interaction in an 
whose track is seen in the central tracker and if the g is also in an state 
either in the inner muon chamber or in the outer muon chamber (it is not 
required that the track be seen in both); b) the momentum of the i?^ track 
in the tracker is measured to be > 10 GeV. c) the energies measured (in an 
appropriate cone surrounding the charged track) by the hadronic calorimeter 
and electromagnetic calorimeter are less than 6 GeV and 2 GeV, respectively 
(both conditions are required to be satisfied, but only the first is relevant for 
a ^-jet). 
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We thank H. Frisch and J. Hauser for clarifying this procedure for us. 
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If an event contains a muonic jet, then the event is discarded in the CDF 
analysis we later employ. Otherwise, the energy of every jet is simply taken 
equal to the energy as measured by the calorimeters. 

• At OPAL|^ the final magnet yoke acts both as the hadron calorimeter and 
the final iron prior to the muon detector. A jet is said to contain a muon if 
there is a charged track in the central tracker, an associated charged track 
in one of the scintillation layers of the hadronic calorimeter and a track in 
the muon chamber. For a ^-jet, we have approximated their procedure by 
requiring that the g be in an state: a) in the tracker; b) as it enters the 
hadronic calorimeter; and c) as it exits the hadronic calorimeter. 

OPAL does not discard events when one or more of the jets contains a muon 
identified in the above way. Rather, the jet energy is corrected assuming that 
the charged track identified as a muon is, indeed, a muon. The procedure for 
computing the jet energy is as follows. 

— Four-momentum vectors are formed for each track and calorimeter clus- 
ter to be included in the jet, and then summed. The three-momentum 
employed for a given track is directly measured in the tracker and the 
energy component for the track is computed by assigning it the pion 
mass, unless it is identified as an electron or muon. (For our purposes, 
we can neglect the e,/x, tt masses.) Calorimeter clusters are treated as 
massless particles; the magnitude of the three-momentum is taken equal 
to the energy of the cluster as measured by the calorimeter. 

— To reduce double counting, four-vectors based on the average expected 
energy deposition in the calorimeter of each charged track are then sub- 
tracted. 

For a ^-jet that has tracks in the tracker and muon chamber that are 
identified as belonging to a muon, this means that the energy and momentum 
vector magnitude assigned to the ^-jet will be given by adding the track 
momentum as measured in the tracker to the total calorimeter response, and 
then subtracting 2 GeV to account for the energy deposit of the supposed 
minimal-ionizing muon. If an R^ track in the tracker does not have an 
associated penetrating track in the muon system (according to the above- 
stated criterion), the track is assumed to be that of a charged pion (it would 
not be identified as an electron), in which case the energy subtracted will be 
taken to be that of a pion with the same momentum as measured for the 
R^ in the tracker. Neglecting the pion mass, this subtraction is equal to the 
measured momentum, with the result that the energy assigned to the ^-jet 
will equal that measured by the calorimeter. Algebraically, we can represent 

^^We thank R. Van Kooten for clarifying the OPAL procedures for us. 
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these alternatives by writing 

Ejet = Pjet = ^Ximeter + ^ (/"id) i^^Pl ' 2 GeV) , (15) 

where 6{fiid) = 1 or according to whether there is or is not, respectively, 
an track identified as a muon associated with the ^-jet. Note that it is 
always presumed that the ^-jet is massless so that Ejet = Pjet is presumed to 
apply. In the OPAL analyses, E^et = Pjet will be defined by this experimental 
procedure and will not be the true jet energy or momentum. 

• A possibly tricky case arises when the R hadron is neutral and undergoes a 
hadronic interaction in the iron of the hadronic calorimeter (or in the unin- 
strumented iron preceding the outer muon chamber at CDF) at a location 
that is less than (roughly) a pion interaction length away from a muon cham- 
ber. This could result in a charged track or, even more probably, a "shower" 
of particles entering the muon chamber from the outer edge of the iron. The 
result would be an anomalous muon signal in the muon chamber. In addi- 
tion, for a track or shower from a hadronic interaction at the edge of the 
hadronic calorimeter, the full energy loss of the i?-hadron from this interac- 
tion would not be measured by the calorimeter. These effects fall outside the 
simplified treatment that we shall employ, described above, which assumes 
that the shower from a hadronic interaction is completely contained in the 
iron. They will be discussed at the end of this section. For now, we present 
results obtained assuming complete containment. 

In order to assess the implications of the OPAL and CDF procedures, we have 
computed the average result for the energy (= momentum), E^^t, assigned to a 
gluino jet for 1000 ^'s produced with a given initial jS, following the OPAL and 
CDF procedures. Since the missing momentum for a given ^-jet is the difference 
between the experimental measurement, Ej^t, and the true initial momentum of 
the g, our focus will be on expectations for the ratio Ejct/ptmc- All results for Ejet, 
here and in future sections, will assume that the shower from a hadronic interaction 
occurring in the iron of the hadronic calorimeter is fully contained. As discussed 
just above, we believe that the effects of incomplete shower containment are small. 

Consider first the CDF detector configuration. We assume A^'cou = 6 interactions 
in instrumented iron and A^coU = 2 uninstrumented interactions between the inner 
muon chamber (which is just outside the hadronic calorimeter) and the outer muon 
chamber. When the gluino is initially produced, and after each subsequent hadronic 
interaction, it is assigned charge \Q\ = 1 with probability P and Q = with 
probability 1 — P. Ionization energy losses are incorporated for any path segment 
between hadronic interactions for which \Q\ = 1. Ionization energy losses are 
multiplied by r = 1.6 when computing the calorimeter response. At each hadronic 
interaction the (AE) of Fig. ^ is assumed to be deposited in the calorimeter and 
included in the calorimeter response (with coefficient 1). If the g is charged in 
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the first track segment, charged after 6 interactions and/or also charged after 8 
interactions, (and has non-zero kinetic energy where it is seen to be charged), and 
the earher described momentum and energy deposit requirements are satisfied, 
then we presume it will be identified as a muon and the ^-jet is discarded. If it is 
not identified as a muon then the ^-jet is retained and the jet energy is set equal 
to the energy as measured by the calorimeter. 




Figure 10: For P = 1/2 and 3/4, we plot, vs. the gluino's initial f3, the average 
fraction of gluino jets that is retained when the CDF procedure is followed. Results 
are given for g masses of m = 5, 25 and 140 GeV, taking r = 1.6. The two curves 
for a given mass are for (AE) cases (1) and (2), the lower curve corresponding to 
case (2). 



The first important issue with regard to the CDF procedure is the fraction of 
^-jets that are discarded as a result of the ^-jet being declared to be "muonic 



)5 



(according to the earlier-stated criteria). In Fig. 10, we plot the average fraction 
of ^-jets retained as a function of the gluino's initial (3, for P = 1/2 and 3/4. 
Results are given for m~ = 5, 25 and 140 GeV. This figure shows that there 
is an intermediate m~-dependent range of (3 for which the ^-jet is "muonic" a 
significant fraction of the time. This occurs as a result of the fact that the energy 
(from electromagnetic and hadronic energy deposits) measured by the hadronic 
calorimeter drops below 6 GeV at intermediate (3. (This happens because, when 
present, the is not sufficiently heavily-ionizing at intermediate j3, and hadronic 
energy deposits typically only become large at large /?.) Note that Fig. [1^ shows 
that events are discarded over a larger range of f3 for (AE) case (2) as compared to 
case (1), in agreement with expectations following from the fact that case (2) yields 
smaller hadronic energy deposits. For P = 0, all ^-jets are, of course, non-muonic 
and are retained. For P = 1/4, the fraction of retained ^-jets is above 0.87 for all 
/? values for all masses and both {AE) cases. P = 1 is a bit of a special 
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we now describe. 

For P = 1, there are no charge fluctuations and, for a given f3 and (AE) case, all 
^-jets are either retained or discarded. For {AE) case (1), we find that the ^-jets are 
retained for all values of (3 for all three m~ values because the hadronic calorimeter 
energy deposits (including both ionization and hadronic collision energy deposits) 
are large enough to fail the < 6 GeV criterion for a muonic jet. For (AE) case (2), 
there is an intermediate range of (3 (dependent upon the value of m~) for which 
the hadronic calorimeter energy deposits are small enough to satisfy the < 6 GeV 
criterion and the ^-jets are discarded as being muonic. These intermediate ranges 
appear as gaps in the (AE) case (2) curves for P = 1 in Fig. ^ below. As a result, 
it turns out that there is a very large difference in the ability of the jets + missing 
energy CDF analysis to exclude a heavy ^-LSP in case (1), which yields good 
sensitivity, as compared to case (2), which yields poor sensitivity. This is clearly 
an artifact of the published CDF analysis procedures. To avoid this sudden change 
in efficiency, we recommend that CDF re-analyze their data without discarding 
muonic jets. 

The second important issue is the measured energy of the retained ^-jets. In 



Fig. |TT|we plot the average (over 1000 produced ^'s) energy assigned to the accepted 
^-jets divided by their actual initial momentum for P = 1/4, 1/2. 3/4 and 1. 
Remarks relevant to borderline cases that will be important in the CDF jets + 
missing momentum analysis are the following. 

• For m~ = 25 GeV and initial j3 > 0.95, the fraction X of the ^'s actual 
momentum that is included in the measured Ej^t is in the range X < 0.15 
for all P values and both (AE) cases. 

• For m~ = 140 GeV and initial (3 > 0.6, one finds X < 0.1 for all P values 
and both (AE) cases. 

The only exception to these generalities occurs when P = 1 and for (AE) case 
(2), for which ^-jets with the above masses and (3 values are discarded as being 
muonic. Aside from this, we can anticipate that gg production at CDF will result 
in an event with large missing momentum. 

In the case of OPAL, if the ^-jet has \Q\ = 1 in the tracker and if it emerges 
into the muon chamber with |Q| = 1 and positive kinetic energy after A^^cou = 7 
interactions then it is assumed that the track in the tracker will be identified as a 
muon and that the jet energy correction of Eq. (|T5|) will be applied. If there is no 
track identified as a muon then the jet energy is set equal to the energy as measured 
by the calorimeter. In Fig. [1^, we plot the average (over 1000 produced ^'s) energy 
assigned to the ^-jet divided by its initial momentum for P = 1/4, 1/2. 3/4 and 
1. For P < 1/2, the f3 ranges of importance at LEP will be those where i^jet is 
only a fraction of the full initial momentum of the g. This is not unlike the CDF 
result. However, for large P there are very substantial differences as compared to 
CDF. For example, when P = 1 most of the kinetic energy is deposited in 
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Figure 11: For P = 1/4, 1/2, 3/4 and 1, we plot, vs. the gluino's initial /3, the 
average measured jet energy E^^t as a fraction of the gluino's initial momentum for 
^-jets that are not declared to be muonic [using the CDF procedures]. Results are 
given for m = 5, 25 and 140 GcV, taking r = 1.6. The two curves for a given mass 
are for (AE) cases (1) and (2). Raggedness in the numerical resuhs, reflecting 
the fact that in our approximation the hadronic interactions only occur at precise 
intervals of 19 cm whereas ionization losses occur continuously, has been smoothed 
out in the plots. Gaps in the case (2), P = 1 curves are where the ^-jet is declared 
to be muonic. 



27 





1.0 




0.8 


Oh 










0.6 




0) 










0.4 












0.2 






0.0 



: OPAL 



P(charged)= 1/4 




:m-=140 



: OPAL 


1 ' ' ' ' 1 ' ' ' ' - 
m = 5 




- ^m — 25 \ 


: m = 140 


1 , , , , 1 , , , , " 



0.4 0.6 0.. 



0.4 0.6 0.. 



P(charged)=3/4 



P(charged) = 1 



Oh 



1.5 



1.0 — 



0.5 



0.0 




Figure 12: For P = 1/4, 1/2, 3/4 and 1, we plot, vs. the gluino's initial (3, the 
average jet energy i^jct [computed using the OPAL procedures, cf. Eq. (^)] as a 
fraction of the gluino's initial momentum. Results are given for m = 5, 25 and 
140 GeV, taking r = 1.6. The two curves for a given mass are for (AE) cases 
(1) and (2). Raggedness in the numerical results, reflecting the fact that in our 
approximation the hadronic interactions only occur at precise intervals of 19 cm 
whereas ionization losses occur continuously, has been smoothed out in the plots. 
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the form of ionization energy losses. If its /3 is too small for penetration to the 
muon detector, then the calorimeter response gives £^jet close to r = 1.6 times the 
g kinetic energy. Once the (3 is large enough for penetration to the muon chamber 
and the tracker track is identified as a muon, £jet, as determined from Eq. ([TSl) , 
jumps to a level that reflects the addition of the g momentum as measured for the 
charged track in the tracker. For P = 3/4 one is in transition from the typical low 
P situation to P = 1. To interpret -Ejet/Ptrue > 1 it is important to recall that it is 
l-^jet ~Ptrue| that determines whether the ^-jet will result in missing momentum. 
Values of E-^ct/ptme significantly different from 1 (whether larger or smaller) will 
lead to missing momentum. Thus, at OPAL, events containing ^'s will generally 
have some missing momentum even when P is large. 

With regard to values of m~ and associated typical /?'s that will be interesting 
borderline cases for the OPAL jets + missing momentum analysis, we note the 
following. 



Consider first m~ = 5 and /? ~ 0.98. Fig. |12] shows that if P is not large, then 
the measured jet energy is small and there will be large missing momentum 
associated with a ^-jet. If P ~ 1, -Ejet/ptrue is somewhat bigger than 1, which 
as noted above will lead to some missing momentum, but not as much as is 
typical at lower P. 



For m~ = 25 and 0.5 < /3 < 0.8, Fig. |I2| shows that the measured jet energy 



is typically a significant fraction of the true momentum once P > 1/2. For 
P = 1, -Ejet/Ptruc is not far from 1 for this (3 range. 

Thus, we can anticipate that P = 1 will yield the weakest OPAL signal at both 
ends of the mass range of interest. 

Hopefully, the discussion of this subsection has provided intuition as to the 
characteristics of ^-jets as measured in the CDF and OPAL detectors. We have 
presented results for what we believe to be the most resonable choice of the inter- 
action length At of the gluino. However, it will be important to assess sensitivity 
to changes in Ay. Smaller At (larger total cross section) yields more hadronic col- 
lisions and, therefore, more hadronic energy deposit and more slowing down of the 
g; larger At, the reverse. We have found that the greatest sensitivity to At arises 
in the case of the CDF jets + missing momentum analysis where larger At implies 
that the smaller hadronic energy deposits and smaller ionization energy deposits 
(due to less rapid slowing down of the g) result in many ^-jets being declared to 
be muonic when P is large, implying a loss of sensitivity for the published analysis 
procedures. In order to provide a representative sample of possibilities for both 
(AE) and At, we will consider three scenarios (denoted SC) in the jets -|- missing 
momentum analyses that follow: 

• SCI: At = 19 cm (as employed in the discussion and graphs given earlier in 
this section) and (AE) case (1). 
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• SC2: At = 9.5 cm and (AE) case (1), implying twice as many hadronic inter- 
actions, and, therefore, larger measured energy for a given ^-jet as compared 
to the SCI case. 

• SC3: \t = 38 cm and (AE) case (2), implying only half as many hadronic 
interactions and small energy deposit per hadronic collision, leading to much 
smaller measured energy for a given ^-jet as compared to the SCI case. 

In the OPAL and CDF analyses of the next sections, our procedure will be to 
generate events containing a pair of gluinos, and then let each gluino propagate 
through the detector allowing for charge changes according to a given choice of the 
probability P at each hadronic interaction. The frequency of hadronic interactions 
is determined by the choice of At, and the amount of energy deposit at each 
interaction is determined by the (AE) case. The characteristics of each event 
are then computed, including overall missing momentum, jet kinematics, etc. The 
relevant cuts are then applied. Only this type of Monte Carlo event-by-event 
procedure allows for all the different types of fluctuations in charge, velocity and 
so forth that take place if gluino-LSP's are being produced. 

3.4 Effects of incompletely contained hadronic interaction 
showers 

Finally, let us now return to the effects that arise if there is a hadronic shower 
at the outer edge of the hadronic calorimeter and, in the case of CDF, at the 
outer edge of the iron shield between the inner and outer muon chamber. This 
mainly affects the jets + missing momentum analyses of OPAL and CDF and 
the heavily-ionizing track analysis of CDF. The details of these analyses will be 
discussed in later sections, but we find it convenient to summarize the influence 
of edge-showers here. We have studied the effects on the analyses in the following 
very extreme approximation. We assume: a) that the last hadronic interaction 
in the calorimeter is completely uncontained and therefore does not contribute to 
measured ^-jet energy; and b) that the last hadronic interaction in the hadronic 
calorimeter, and, for CDF, also the last interaction in the iron shield, yields a 
charged track in the subsequent muon chamber. We flnd the following results. 

• Small P: In the OPAL and CDF jets -|- missing momentum analyses, the 
jet is declared to contain a muon only if a charged track is also seen in the 
tracker. For small P, this probability is small. The main effect would then 
be that the energy of the hadronic interaction shower at the edge of the 
calorimeter would not be deposited in the calorimeter, thereby leading to a 
decrease in the measured jet energy. We find that the resulting increase in 
missing momentum would be modest 10 — 15%), even in our extreme 
approximation. This would yield some enhancement in the efficiency for the 
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jets + missing momentum signal in the OPAL and CDF analyses, but not 
enough to significantly alter the limits on m~ that are obtained. 

The heavily-ionizing track signature is not relevant for small P since there is 
low probability for a charged track in the tracker. 

• Large P: For large P values, in the jets + missing momentum OPAL analy- 
sis, the ^-jet will be declared to contain a muon regardless of whether there 
is an extra muon-chamber track or shower. Also, since most of the i?-hadron 
energy losses are in the form of ionization rather than from hadronic interac- 
tions, we find that the measured ^-jet energy only decreases slightly. Thus, 
the OPAL jets + missing momentum results would be little affected. 

Turning to the CDF jets + missing momentum analysis, we again note that, 
when P is large, most of the measured energy is from ionization energy 
deposits and earlier hadronic interactions, and the incomplete containment of 
the tracks/shower of a last hadronic interaction in the hadronic calorimeter 
generally has little affect, provided the ^-jet is declared not to be muonic. 
(Note that if the incompletely contained shower originates in the outer edge 
of the iron between the inner and outer muon chambers it would not have 
been instrumented, i.e. would not contribute to measured energy anyway.) 
Unless one is right on a borderline, the small decrease in measured energy 
due to losing the shower from the last hadronic interaction in the calorimeter 
will not cause a ^-jet that would otherwise be declared to be non-muonic to 
fall into the muonic category. However, we have already seen in Fig. |TT] that 
for P = 1 we are right on such a borderline, with case (2) {AE) giving rise 
to large gaps (in (3) for which the ^-jet is declared to be muonic whereas for 
our SCI case (1) choice the ^-jet is never declared to be muonic. We find 
that failure to capture any of the energy of the last shower also pushes us 
past this borderline. Thus, in our extreme approximation, the loss of the 
shower results in much the same phenomenology for CDF as the SC3 case 
defined earlier; one finds that a very substantial weakening of the jets -|- 
missing energy signal occurs. Of course, as already noted earlier, the way 
around this is to re-analyze the CDF data without throwing away muonic 
jets, perhaps using something like the OPAL procedure. 

• Moderate P: For moderate P values, the penetration of a hadronic interac- 
tion shower to the muon chamber would tend to increase the number of ^-jets 
that are declared to contain a muon in the OPAL analysis. The momentum 
computed for the extra muon-jets via Eq. ([151) will be substantially larger 
than otherwise. On average this increase in momentum is only partially off- 
set by the decrease in the measured calorimeter energy deposit from the jet 
due to non-containment of the final shower in the hadronic calorimeter. The 
net result is a modest decrease in the efficiency for the jets -|- missing energy 
signal. However, the m~ limit borderline is so sharp at moderate P (see later 
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OPAL results) that there would be little change in the limits that can be 
extracted from the OPAL analysis. 

In the CDF analysis, there are two effects. The extra muon-chamber signal 
will tend to decrease the number of non-muonic events because a) there are 
more events with tracks in the muon chambers and b) because the energy 
deposit measured by the hadronic calorimeter decreases as a result of in- 
complete containment of the tracks of the final shower. However, a sizeable 
fraction (roughly, 50% for (AE) case (1) and P = 1/4, 1/2, and 3/4, in the 
f3 regions of relevance) of the events that are retained at moderate P (see 



Fig. |IOD are non-muonic because of the absence of a charged track in the 
tracker. The retention of these events would be unaffected by the presence of 
an anomalous muon-chamber signal. Overall, we find that the decrease in the 
number of accepted ^-jets is typically of order 30%. However, this decrease 
is compensated by the fact that the decrease in measured calorimeter energy 
due to incomplete shower containment increases the missing momentum and, 
therefore, the efficiency for non-muonic events that contain such a shower. 
(Recall that, once accepted, the ^-jet momenta are computed in the CDF 
analysis without including any muon correction.) Changes in the extracted 
m~ limits would not be large. 

• For moderate or large P: The heavily-ionizing track (HIT) searches that can 
be used to eliminate a span of m~ values when P > 1/2 will be completely 
unaffected by an anomalous muon-chamber signal in the case of OPAL (since 
the OPAL HIT analysis, described later in section 6, essentially only uses 
tracker information) and will be enhanced in the case of CDF (since the CDF 
HIT analysis, discussed in section 7, requires a track in the inner and/or outer 
muon chamber in addition to a HIT in the inner tracker). 

Thus, we think that the effects upon our analyses of a hadronic collision that 
leads to an anomalous muon-chamber track or shower are small, except in the 
case of large P in the jets -|- missing momentum CDF analysis where one is very 
sensitive to just how much of the energy in the final hadronic calorimeter shower 
escapes into the muon chamber. We repeat our expectation that this sensitivity 
could be eliminated by removing the "non-muonic" jet requirement in the CDF 
analysis. A study of the effects of incomplete shower containment is probably best 
left to the detector groups themselves. 

Finally, we note that events having a shower entering the muon chamber would 
actually appear to provide a potentially spectacular signal for a ^-LSP — one 
that should be specifically searched for. This signal would appear to be especially 
promising if P is small and one focuses on events in which there is no charged track 
in the tracker associated with the jet pointing to the muon chamber shower. 
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4 Constraints from LEP and LEP2 



At LEP and LEP2, we assume that all other SUSY particles are beyond the kine- 
matic reach of the machine. The only possible signal for SUSY is then pair produc- 
tion of two gluinos. Gluinos can only be produced via two processes: e'^e'^ 1^99 
271 PSj p9| , which can take place at tree- level, and e"'"e~ — > gg 0, which 



takes place via loop diagrams (involving squarks and quarks). As discussed later, 
the latter process is very model dependent and can be highly suppressed. Thus, we 
begin by focusing on the qqgg final state. We consider both the LEP Z-pole data 
and higher energy running at LEP2. The (uncut) qqgg cross section^ is plotted 



in Fig. |T3| as a function of m~ for ^/s = mz, 172 GeV, 183 GeV and 192 GeV. 



Given that the total e+e^ —>■ Z cross section is ~ 6 x 10^ pb. Fig. |T3| implies that 
BR{Z — > qqgg) > few x 10^^ for m~ < 25 GeV. Since millions of Z's have been 
produced at LEP, we can demonstrate that ^'s lighter than this and heavier than 
about 5 GeV can be ruled out. In contrast. Fig. makes it clear that very sub- 
stantial luminosity at higher LEP2 energies will be required for constraints from 
LEP2 data to be competitive. For example, L = 500 pb~^ at y/s = 192 GeV 
will yield only about 4 e+e~ qqgg events (before cuts) at m~ = 25 GeV. Also 



shown in Fig. |T3| is the uncut e~^e —>■ qqgg cross section at ^/s = 500 GeV, a 



possible choice for the next linear collider (NLC). One finds ci{qqgg) < 1 fb for 
m~ > 60 GeV, which would correspond to 50 events for L = 50 fb~^. Even for 
L = 500 fb~^ one finds fewer than 5 events [c^qqgg] < .01 fb] for m~ > 140 GeV, 
which will turn out to be close to the lower limit that can already be set by using 
Tevatron data. 

Thus, we focus on ^/s = mz- The procedures for employing LEP Z-pole data 
to place constraints on the ^-LSP scenario depend upon the manner in which 
the ^-jet is manifested in the detector; this was outlined in the previous section. 
Generally speaking, qqgg events will have 4 jets and missing momentum. As noted 
in the previous section, the most crucial kinematical aspect of the ^-jets is their 
distribution as a function of (3. The number of ^-jets as a function of (3 is presented 
in Fig. for m~ = 5 GeV and 25 GeV. We see that a light gluino with m~ < 
5 GeV has a (3 distribution that peaks at /? ~ 0.98 while a heavier gluino with 
m~ ~ 25 GeV has a broad /? peak centered about /? ~ 0.6, with the most probable 
/3 values lying between 0.5 and 0.7. The implications of these /? ranges at these two 
masses were already indicated in the previous section. The reason that we will not 
be able to obtain limits from LEP data for very small m~ values is that as the gluino 
bound state mass decreases below 5 GeV, the initial (3 of the g increases. As a 
result, the energy loss in the first few hadronic collisions increases significantly. For 
a mass of < 1 GeV, the energy loss is essentially complete (that is the calorimeters 



^^We have employed a numerical helicity amplitude computation for e+e — > qqgg valid for 
itrary m~; the program is availal 

element can also be found in Ref. |32 



arbitrary m~; the program is available upon request. A crossed version of the squared matrix 
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Figure 13: cr(e+e qqgg) as a function of m~ for y/s = mz (solid), 172 GeV 
(dashes), 183 GeV (dot-dash), 192 GeV (dots) and 500 GeV. No cuts. 



will contain the hadron). 

The most relevant LEP experimental analyses currently available are those 
related to the search for pair production of neutralinos, Z —>■ X1X2? with X2 
The OPAL Q and L3 || 



analyses have the highest statistics and place limits 
on X1X2 production in the jets +j/t channel that are potentially relevant for the 
qqgg final state. However, the L3 analysis is restricted entirely to 2j +j/t final 
states. Only the OPAL analysis is relevant to any nj +'^t final state with n > 2. 
Typically, qqgg events give n = 2, 3, or 4, depending upon the amount of energy 
deposition by the ^-jets. 

The OPAL analysis is based on dividing the event into two hemispheres as 
defined by the thrust direction of the visible jets. We have implemented their 
procedures in a parton-level Monte Carlo and computed the efficiency for the Z 
qqgg events to pass their cuts as a function of m~ for various choices of the charged 
fragmentation probability P. Our precise procedures are as follows. In the OPAL 
analysis of multi-jet events, each event is divided into two hemispheres by the plane 
normal to the thrust axis, where the thrust T is defined as 



T = max^ 



Ei \Pi ■ n\ 

Ei \Pi\ 



(16) 



and the thrust axis is the h that leads to the maximum. In the OPAL analysis. 
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Figure 14: Distributions of the number of s^-jets as a function of (3 at LEP {^/s — 
mz) for m~ — 5 and 25 GeV. No cuts are imposed. 

the Pi are assigned to calorimeter clusters and associated tracks as described in 
the previous section. Associated energies are computed as if the track/cluster 
composites have very small mass. The sum of the (visible) four-momenta in a given 
hemisphere defines the four- momentum of the 'jet' associated with that hemisphere; 
note that the 'jet' need not have zero invariant mass. OPAL then separates events 
into mono- or di-'jet' events, where a mono-'jet' event is one having a 'jet' in only 
one hemisphere. Mono-'jet' events are discarded. The following cuts are then 
applied to the di-'jet' events: 

\{M^^- ' + M^^- 2) < 20 GeV , M,,JE,^ > 0.27 , 
Pt > 10 GeV , < 20 GeV , 
T > 0.7 , min[Thcm.i, rhcm.2] > 0.7 , 
cos6'acoi < 0.98 , I cos6'miss| < 0.94 , 
cos^acoi < 0.95 , cos^acop < 0.98 if both 'jet's are in | cos^^| < 0.71 , 
cos^acoi < 0.90 ,cos^acop < 0.95 if either 'jet' is in | cos^| > 0.71 , 

where (tt — ^acoi) is the three-dimensional angle between the two 'jet's, (tt — ^acop) 
is the angle between the two 'jet's in the x-y plane, ^miss is the polar angle of the 
missing momentum, M^is is the visible mass, and p (used to compute p^ and p^) is 
the vector sum of all (visible) three momenta. In the above, M^jg is computed by 
summing all the visible four-momenta (as defined earlier) in the event and taking 
the square. The square of M^^- for each hemisphere is computed by summing 
the visible four-momenta in the hemisphere and squaring. The thrust. Them., for 
each hemisphere is defined by going to the center-of-mass for that hemisphere 
(defined by the sum of all visible three-momenta in the hemisphere being zero) 
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and computing the thrust as in Eq. (|T6D using only the three-momenta of that 
hemisphere. 

In applying the above procedures to the Monte Carlo events, it is necessary to 
adopt an algorithm for including the effects of detector resolution. In our com- 
putations, all cluster/track momenta and energies are smeared using the stated 



OPAL hadronic calorimeter energy resolution of / E = 12Q%/ ^ E[ GeV). We 
note that energy smearing is important in that it generally increases the OPAL 
acceptance efficiencies by virtue of the fact that, on average, jet-energy mismea- 
surement tends to enhance the amount of missing momentum. This enhancement 
is especially important for m~ and P choices {e.g. m~ = 25 GeV and P = 1) such 
that the missing momentum before smearing is small. Another important ingredi- 
ent is properly accounting for the fact that the i?-hadron does not take the entire 
momentum of the g. We have employed the standard Peterson form for the 
fragmentation function oig^R: 



9 



1 1 '5 



z 1 — z 



(17) 



where we will take e~ = (0.3 GeV/m~)^. Here, the i?-hadron carries a fraction 
z of the momentum of the g and a normal (light quark or gluon) jet carries the 
remainder. The i?-hadron is then treated in the calorimeter as we have described 
in the previous section. The energy of the remainder (effectively zero-mass) jet 
is taken equal to its momentum and is assumed to be entirely deposited in the 
calorimeter. Typically, fragmentation does not have a large influence on the effi- 
ciency with which events are retained, especially in cases for which the ^-jet energy 
is measured to be a large fraction of the ^'s initial kinetic energy. 

The OPAL data corresponds to A^had = 4.4 x 10® hadronic Z decays. The 
expected number of qqgg events after cuts is then 

Y - ^hi,dBR{Z qqgg) x efficiency ^^^^ 
BR[Z —>■ hadrons) ' 

where we use the efficiency as computed via the Monte Carlo. After cuts, OPAL 
observes 2 events with an expected background of -B = 2.3 events. The 95% upper 
limit on a possible new physics signal is then S = 4 events, corresponding to 
BR{Z —>■ qqgg) x efficiency ~ 6.4 x 10~^. How low a value of m~ can be eliminated 
depends upon the efficiency at low m~. Because of the very high raw event rate 
at low m~ values, quite small efficiency can be tolerated. We will see that we can 
exclude gluino masses above 3 — 4 GeV. 

As described in the previous section, to obtain a reliable result for the range 
of m~ that the OPAL analysis excludes, we have computed the efficiency for qqgg 
events to pass the full set of cuts when Eq. ([TsD is employed for each g on an 
event-by-event basis, including (for P 7^ 0, 1) random changes (with probability 
determined by P) of the i?-hadron charge at each of the hadronic interactions it 
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Figure 15: In the upper window, we plot the OPAL qqgg event efficiency (after 
all cuts) in the P = 0,1/4,1/2,3/4,1 cases, as computed using event-by-event 
determination of Ej^t [using Eq. (|15D] for each g. For P 7^ 0, 1, changes of the R- 
hadron charge as it passes through the detector are randomly implemented. Both 
smearing and fragmentation effects are included. The lower window gives, as a 
function of m~, the corresponding 95% CL upper limits compared to the theoretical 
prediction for BR{Z qqgd)- Results are for the SCI choices of At = 19 cm and 
(AE) case (1). 
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experiences as it passes through the detector. We have considered the three sce- 
narios — SCI, SC2, and SC3 — for choices of and the (AE) case that were 
outhned at the end of the previous section. In Fig. we plot the resuhing OPAL 
efficiency for qqgg events after all cuts as a function of m~ for P = 0, 1/4, 1/2, 3/4, 1 
for the SCI choices, including calorimeter energy smearing and fragmentation ef- 
fects. Also shown are the resulting 95% CL upper limits on BR{Z qqgg)- We 
see that for any P not near 1, the entire range from low m~ ~ 3 GeV to high 
m~ ~ 25 GeV is unambiguously excluded. For P ~ 1, the largest value of m~ that 
can be excluded is about 23 GeV. [The m~ > 23 GeV limit for P = 1 is similar 
to, but somewhat higher than, the limit obtained by searching for heavily-ionizing 
tracks at OPAL (discussed later in section 6).] 
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Figure 16: 95% CL upper limits as in Fig. |T5| except that we do not include the 
effects of energy smearing or fragmentation. 



In Fig. |T6| we present the 95% CL limits obtained without including either 
energy smearing or Peterson fragmentation. This figure shows that the limits are 
little altered except for P ~ 1, in which case the OPAL analysis does not exclude 
any significant range of m~. It is energy smearing that is the dominant factor in 
obtaining a significant efficiency for event acceptance when P ~ 1. Even though 
P ~ 1 leads to Ej^t ~ Ptrue at the parton level [for the P values typical for the 
m~ =5 — 25 GeV mass range (see Fig. 0)] and thus small missing momentum 
at the parton level, energy smearing produces large event-by-event fiuctuations in 
the measured energy of each g jet which lead to substantial missing momentum for 
many events. 
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Figure 17: 95% CL upper limits as in Fig. 15 except that we use the SC2 choices 
of At = 9.5 cm and (AE) case (1) in the upper window and the SC3 choices of 
At = 38 cm and (AE) case (2) in the lower window. 
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Results analogous to those obtained for the SCI choices of Ar = 19 cm and 
(AE) case (1), and presented in Fig. |15|, are presented for the SC2 and SC3 choices 
[SC2: At = 9.5 cm, (AE) case (1). SC3: At = 38 cm, {AE) case (2)] in Fig. |T7. 
In fact, these possible extremes always give higher efficiencies and a slightly larger 
range of m~ exclusion than found in the SCI case. 

We expect that re-analysis of the LEP data sets using cuts more appropriate to 
the qTjgg final state for given values of P and m~ will yield only a small improvement 
over the results obtained using the existing X1X2 analysis cuts. At large m~, the 
event rates are falling so rapidly that the 95% CL upper limit is not likely to 
be increased by more than a few GeV. Ruling out m~ values significantly below 
3 — 5 GeV will be difficult since for such m~ the gluino looks so much like a normal 
jet that only the still controversial analyses of Ref. [0] are likely to prove relevant. 
Still, we would recommend attempting to make use of the threshold in the mass 
recoiling against the two energetic jets of the the event present at Mrecoii ~ 2m~. 
Perhaps the background could be reduced to zero by an appropriate set of cuts 
including one requiring Mrccoii ^ 2m~. 

It is also worth nothing that the jet energy as computed using the OPAL 
procedure of Eq. ([T5|) is often larger than the actual g energy for large P. This 
may be interesting at LEP, since there it is possible to compare the total measured 
or 'visible' energy associated with an event to the total center of mass energy. By 
summing the assigned energies of all jets, one would find events in which the total 
energy exceeds the center of mass energy when P is near 1. Indeed, the above 
Monte Carlo generates a significant number of such events when m~ is small. To 
our knowledge, the LEP experimental groups have not analyzed their events in a 
manner that would be sensitive to such a discrepancy. 

Finally, we briefiy discuss e~^e~ gg production via quark-squark loops. 
Again, only the existing Z-pole data might possibly yield a useful constraint. As 



discussed in Refs. even if the squarks are all completely degenerate, the 

Z —>■ gg branching ratio can be non-zero by virtue of the top mass being much 
greater than the bottom mass. However, Ref. |^| finds BR{Z gg) < 2 x 10~^ 
for all m~ if the common squark mass (m) is above ~ 200 GeV. The typical event 
would contain two back-to-back jets. But these would not generally have equal 
energy due to the fact that fiuctuations would be substantial, especially if P is in 
a range such that there would sometimes, and sometimes not, be a charged track 
identified as a muon contained in one or both of the jets. For small deposited 
energy per ^-jet, as typical for small P, the net apparent energy of the typical 
event would be below mz, possibly causing such events to be confused with the 
two-photon background. For large enough P and smaller m~, many of the events 
would be anomalous in that the sum of their apparent energies would exceed mz- 
We are uncertain if any of the LEP analyses would have been sensitive to such 
events appearing at a level corresponding to BR{Z gg) ~ 1 — 2 x 10~^. In 
any case, the gg event rate can be suppressed to an unobservable level simply 
by taking m sufficiently large. (Roughly, BR{Z gg) falls as 1/m^.) Thus, no 
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model-independent m~ limits from the gg final state are possible. 

5 Present and Future Tevatron Constraints from 
jets + 

In the ^-LSP scenario, with all other SUSY particles taken to be much heavier, the 
only standard hadron-coUider SUSY signal is jets+j/T- Current mSUGRA analyses 
of this channel do not apply since the g does not cascade decay {g qqXii ■ ■ ■) to 
additional jets. In the ^-LSP scenario, for a given value of m~, fewer hard jets are 
expected and the amount of missing momentum is typically smaller. Consequently, 
the limits that can be placed on m~ from Tevatron data will be weaker. Still, we 
will find that substantial constraints can be placed on the ^-LSP scenario using 
existing Tevatron data, and that even stronger constraints will arise from Runll 
data. 

In assessing the ability of the Tevatron to discover or exclude a heavy ^-LSP, we 
have employed cuts that mimic those employed by CDF in analyzing RunI data in 
the jets+^^T channel. CDF cuts [^, ^] are employed rather than DO cuts |^ since 
the CDF jet-energy andj/y requirements are weaker than required in the DO cuts. 
For the same integrated luminosity, weak cuts allow greater sensitivity to the heavy 
^-LSP situation in which the most energetic jets come from gluons radiated from 
the initial state colhding partons. The precise CDF cuts used are those employed 
in Ref . [|] ; they are designed to duplicate the experimental procedures of Ref . |^ 
to the extent possible in the context of a Monte Carlo simulation. 

• LI: No (isolated) leptons with Et > 10 GeV. 

• MPT:]/t > 60 GeV. 

• NJ: There are n{jets) > 3 with |?7jet| < 2 and Et > 15 GeV, using a coales- 
cence cone size of AR = 0.5. 



Azimuthal separation requirements as follows: 

- JIMPT: A0fe/T,ji) < 160°; 

- JMPT: A(f)^T,j{ET > 20 GeV)) > 30°. 

These are designed, in particular, to reduce QCD jet mismeasurement back- 
ground. 



-'^^The situation being considered is not dissimilar to the O-II model case where the gluino, 
Xi and xt nearly degenerate with one another. The RunI Tevatron limits for this latter 

scenario were determined in Ref. n. 
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Events were generated using ISAJET-7.37 ||3^. Each event was passed through a 
toy calorimeter with cells of size Arj x A0 = 0.1 x 0.1 extending out to \ri\ = 4. 
Electromagnetic and hadronic resolutions of 15%/ and 70%/ V^, respectively, 
were chosen to approximate those of CDF. The most important cut is the missing 
transverse momentum (MPT) cut. This is especially true at low m~. Typically 
only a small fraction of the events are retained after the MPT cut. The next most 
important cut is the jet-number (NJ) cut. Typically, for P and m~ choices that 
give larger MPT cut acceptance, the NJ cut acceptance is smaller. At the higher 
140 GeV mass, the cuts retain a larger fraction of events than at lower mass. 



m 



(But, of course, the cross section is smaller at high mass.) 
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Figure 18: The (3 distributions of the g^s produced in pp gg, before cuts, for 
m~ = 40 GeV and m~ = 140 GeV, taking = 1.8 TeV. 

In order to relate the Tevatron situation to the discussion of section 3, it is 
useful to present the (3 distribution of the g for several m~ values. In Fig. [1^, we 



present the (3 distributions, before cuts, for m~ = 40 GeV and 140 GeV, i.e. values 
near the upper and lower ends of the interesting mass range. For m~ < 40, /3 is 
typically > 0.95; for m~ ~ 140 GeV, the (3 distribution peaks near /3 ~ 0.75, with 
most events having 0.5 < [3 < 0.9. The /5 distributions, both before and after cuts 
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Figure 19: In the top window, the (3 distribution of the gf's produced in pp — > gg, 
before cuts, for m~ = 100 GeV, taking ^/s = 1.8 TeV. In the lower three windows, 
distributions in /3 after cuts are compared for P = 0, 1/2 and 1. 
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(taking P = 0, 1/2 and 1), are given for m~ = 100 GeV in Fig. |T9|. Referring 
back to Fig. |ll| and related comments, we see that in all cases the most probable 
(3 values are such that the measured E^^t of most ^-jets will be much smaller than 
the true momentum, thereby leading to large missing momentum as defined in the 
analysis. 
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Figure 20: Thej/x distribution (before cuts) for pp gg events at ^/s = 1.8 TeV 
is illustrated for m~ = 100 GeV and P = 0, 1/2 and 1. 



The distribution in^T that results is illustrated for m~ = 100 GeV in Fig. pO 



There, we see a substantial tail with^j- > 60 GeV that is essentially independent 
of the choice of P. This independence of P is due to the small dependence of the 
/3 distribution on P (as illustrated in Fig. ITW) and to the CDF procedure in which 
events where one of the ^-jets looks muonic are discarded and no correction is 
applied to the calorimetric energy measurement for a retained ^-jet that contains 
a penetrating track. 

Let us now turn to determining the limits on a ^-LSP from the CDF data. 
To do so, we compare the cross section for gg pair production after cuts to the 
SM background expected by CDF. For the above CDF cuts and y/s = 1.8 TeV, 
Ref. quotes a background rate of 28.7 events for L = 19 pb~^, corresponding 
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to (Tb = 1-51 pb. (A background rate of 33 events is quoted for the very slightly 



different > 3 jet cuts of the final published CDF analysis, Ref. we prefer to 



stick to the cuts of Ref. [^].) The 95% CL lower limit on m~ is obtained when 
the signal rate declines below the 1.96cr level, corresponding to as ~ 553 fb (after 
cuts). We note that this is about the same as the as ~ 614 fb required for a 5a 
signal at L = 0.1 fb~^. This latter cross section level will be indicated on our 
figures. In Runll, systematic uncertainties in the background will very probably 
determine the limit of sensitivity. Indeed, the 95% CL and 5a levels for as are 
much lower for L >2 fb^ than the as sensitivity limit defined hy S/B > 0.2 {i.e. 
as > 302 fb). For instance, the 95% CL cross section upper limits would be 53.9 fb 
(15.2 fb) for L = 2 fb~^ (25 fb~^), respectively. If systematics can be understood 
at a better than 20% level, then the limits that could be obtained from Runll 
using RunI cuts would improve substantially as compared to the S/B > 0.2 level 
limits. Correspondingly, a ^-LSP signal with S/B 0.2 would have a very high 
nominal S/ \fB. Clearly, optimization of the cuts and procedures can be expected 
to improve upon these first estimates of sensitivity at Runll. 



In Figs. ^ |2^ and|2^, we plot the cross section, as, after cuts, as a function 
of m~ for P = 0, 1/2, 3/4 and 1, for the SCI choices of At = 19 cm and {^E) 
case (1). Also shown on these plots is the L = 0.1 fb^^ S/\fB = 5 cross section 
level (which, as discussed above, is about the same as the 95% CL lower limit for 
L = 0.19 fb-i). We see that, at 95% CL, current CDF analyses [||, of the 
L = 19 fb"^ data set require m~ > 150, 130, 130, 140 GeV for P = 0, 1/2, 3/4, 1, 
respectively, and that, for all P, m~ values are excluded from the upper limit all the 
way down to < 20 GeV at a very high CL. Note that the 130 — 150 GeV lower limit 
on m~ obtained is substantially below the lower limit that RunI data places on m~ 
in a typical MSSM model. For easy comparison. Figs. |2^ and |2^ all show 

the cross section (after cuts) resulting from gluino pair production in the MSSM 
model considered in Ref. [^] with m~ = 1000 GeV, /x = —400 GeV and tan /? = 4; 
one sees that RunI data yields a 95% CL limit of roughly m~ > 210 GeV. 

We re-emphasize that in the Monte Carlo we have treated each ^-jet on an event- 
by-event basis. In this way, the decision as to whether a given ^-jet is "muonic" 
is made event-by-event, including (for P < 1) the possibility of charge changes 
(allowed for in random fashion on an event-by-event basis according to the chosen 
P) at each hadronic interaction as the g traverses the detector. 

As for the OPAL analysis, we wish to assess sensitivity of our CDF results to 
the choices of Xt and (AE) case. In order to do so we present several results for 
the extreme choices defined earlier in section 3, and denoted by scenario labels SC2 
and SC3. First, in Figs. pB|and we present P = 3/4 results for the SC2 and SC3 
choices, respectively. We observe that when P is large SC2 (SC3) choices result 
in stronger (much weaker) limits from the CDF analysis. The poor SC3 results 
are easily understood as follows. For the SC3 choices, significantly less energy is 
deposited by a ^-jet. (The hadronic energy losses are smaller for the longer At and 
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Figure 21: The cross section (after cuts) in the jets+^^r channel is compared to (a) 
the 5(7 level for L — 0.1 fb""^ (also roughly the 95% CL upper limit for L — 19 pb~^) 
at = 1-8 TeV and (b) the S/B = 0.2 level at Runll (L > 2 fb"\ = 2 TeV) 
as a function of m~ for P — 0. SCI choices of At = 19 cm and (A£') case (1) are 
employed. 
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Figure 22: The cross section (after cuts) in the jets +j/t channel is compared 
to (a) the 5(T level for L = 0.1 fb^^ (also roughly the 95% CL upper limit for 
L = 19 pb"^) at ^/s^ 1.8 TeV and (b) the S/B = 0.2 level at Runll (L > 2 fb-\ 
a/s — 2 TeV) as a function of m~ lor P — 1/2, using event-by-event determination 
of the momentum (=energy) of each ^-jet (including the probabilistic treatment 
of charge-exchanges at each hadronic collision) in events such that neither ^-jet is 
"muonic" (see text). SCI choices of At = 19 cm and {^E) case (1) are employed. 
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Figure 23: The cross section (after cuts) in the jets +j/t channel is compared 
to (a) the 5(T level for L = 0.1 fb^^ (also roughly the 95% CL upper limit for 
L = 19 pb"^) at ^/s^ 1.8 TeV and (b) the S/B = 0.2 level at Runll (L > 2 fb-\ 
a/s — 2 TeV) as a function of m~ for P = 3/4, using event-by-event determination 
of the momentum (=energy) of each ^-jet (including the probabilistic treatment 
of charge-exchanges at each hadronic collision) in events such that neither ^-jet is 
"muonic" (see text). SCI choices of At = 19 cm and {AE) case (1) are employed. 
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Figure 24: The cross section (after cuts) in the jets+^/r channel is compared to (a) 
the 5(J level for L = 0.1 fb~^ (also roughly the 95% CL upper limit for L = 19 pb~^) 
at = 1.8 TeV and (b) the S/B = 0.2 level at Runll (L > 2 fb"\ ^ = 2 TeV) as 
a function of m~ for P = 1, using event-by-event determination of the momentum 
(=energy) of each ^-jet in events such that neither ^-jet is "muonic" (see text). 
SCI choices of At = 19 cm and (A£') case (1) are employed. 
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Figure 25: As in Fig. 23, except that SC2 choices of Xt = 9.5 cm and {AE) case 
(1) are employed. 
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Figure 26: As in Fig. 23, except that SC3 choices of Xt = 38 cm and {AE) case 
(2) are employed. 



51 



o 
b 



o 
u 

0) 

-p 
b 



10^ 



10^ 



10^ 



10' 



10 



lo-- 



10^ 



10^ 



10' 



10- 



FT \ I \ ^ I \ I \ I \ I \ \ ^ \ \ \ \ ^ \ I \ r 

V^=1.8 TeV 

t ' + * 

+ X 

X MSSM 



- 5o" limit for 0.1 fb' 



g=LSP + 
P=l; SC2 + 



J \ I \ \ I \ L 



50 100 150 200 250 



- 1 1 1 1 1 


1 1 1 1 1 


1 1 1 1 1 1 1 1 z 

^=2 TeV : 




* 

+ 


X I 


g = 


LSP 


X MSSM - 


p= 


1; SC2 


X 

+ X 






+ = 

X 


- S/B>0.2 


for 2 and 


25 fb"^ + 

+ - 

1 1 1 1 1 1 1 1 


1 1 1 1 1 


1 1 1 1 1 



50 100 150 200 250 

(GeV) 



Figure 27: As in Fig. 24, except that SC2 choices of Xt = 9.5 cm and {AE) case 
(1) are employed. 
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Figure 28: As in Fig. 22, except that SC3 choices of Xt = 38 cm and {AE) case 
(2) are employed. 
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smaller case (2) {AE)^s, and the ionization energy losses are smaller because the 
g does not slow down as much due to the smaller hadronic energy losses.) As a 
result, when P is large the ^-jet is much more likely to be declared to be "muonic" , 
both because it is highly probable that it will make it to either the inner or outer 
muon chamber, and be charged therein, and also because the total energy deposit 
will not exceed the CDF cutoff and thereby prevent its being declared to be a 
"muonic" jet. Thus, many more events are discarded. As P increases above 3/4, 
the cross section obtained for a given m~ after cuts decreases further. For example, 
for P = 1 current CDF jets + missing momentum data and analysis procedures 
provide no constraints on m~ for the SC3 choices, whereas Fig. ^ shows that 
strong constraints are provided for the SC2 choices. Finally, in Fig. we show 
that, for P = 1/2 (and smaller), even if we make the SC3 choices the limits on 



m~ are nearly as strong as for the SCI choices of Fig. 22. For SC2 choices, the 



corresponding plot would show even stronger limits than for the SCI choices. 

Thus, the jets + missing momentum data and analysis of CDF only allows a 
g with m~ < 130 GeV if the g has a high charged-fragmentation probability and 
rather weak hadronic interactions. Fortunately, the CDF heavily-ionizing-track 
analysis discussed later provides strong constraints for large P that exclude this 
possibility for m~ > 50 GeV (which should be extendable to lower m~ values). As 
we have repeatedly noted, the lack of sensitivity of the RunI CDF jets + missing 
momentum analysis would disappear if the data is re-analyzed without eliminating 
events containing a muonic jet. We urge the CDF collaboration to perform this 
re-analysis. 



As one possible backup at low m~, we looked at whether or not UAl |^ and 
UA2 |5T] data could be used to exclude m~ in the m~ ~ 30 GeV region. We find, 
however, that no limits on m~ in this (or any other mass region) are possible from 
the UAl and UA2 data. Another backup at low m~ could be an analysis of pre- 
scaled data {i.e. data not taken at the full trigger rate) accumulated using lower 
Pt cuts on the jets. For example, CDF took about 1 pb^^ of data using a \ow-Et 
four-jet trigger Such data might be useful since at lower m~ the standard 

CDF cuts employed above tend to yield a rather small efficiency for accepting 
signal events. We have not examined this data in detail. 

Let us now consider Runll. Returning to Figs. 0, ^ and we see that 
the limits based on S/B > 0.2 will rise to m~ > 180, 160, 160, 180 GeV for P = 
0, 1/2,3/4, 1, respectively, for Runll (with L > 0.5 fb^^). If systematics could be 
controlled so that a signal with S/B < 10% becomes reliable, each of these lower 
limits would be increased by about 30 GeV. All these potential lower bounds are, of 
course, still substantially lower than the m~ lower bound that can be achieved in the 
reference MSSM model for the same S/B criterion {e.g. 250 GeV for S/B > 0.2). 
It is worth noting that Runll limits will be much less sensitive to At and (AE). As 
shown in Fig. even the SC3 choices will allow exclusion of all m~ < 130 GeV. 

We end by noting that if the squarks are not much heavier than the g, then the 
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gg cross section at the Tevatron will be reduced due to negative interference effects 
in the qq — gg amplitude from squark exchanges. However, the gg gg amplitude 
is unaffected. Further, additional very prominent signals will emerge from squark 
production channels that will more than compensate. Thus, the approach of taking 
all other SUSY particles to be much heavier than the g can be expected to yield 
the most conservative limits for the ^-LSP models. 

6 The OPAL signal for a charged gluino hadron 



OPAL has searched for e+e^ QQgg events in which the ^'s fragment to a 
charged that traverses their two-meter radius tracking chamber. They look 
for events with an anomalous value for the ionization dE/dx as compared to the 
momentum Both quantities are measured in the tracking chamber. As a result, 
penetration of the track to the muon detectors is not required. After appropriate 
kinematical cuts and cuts on the region of the dE/dx — |p| plane accepted, there 
is only one candidate event. They convert this into a 95% CL limit on the number 
of signal events. To interpret this limit they compute the expected number of 
gluinos produced and accepted and multiply by the probability P for g — > 
fragment at ion. 13 They place 95% CL upper limits on P as given in Table |l[ 



m~ 


1.5 2.3 3.0 5.0 10.0 15.0 20.0 


p max 
-'^95% CL 


0.37 0.20 0.14 0.06 0.13 0.33 1.03 



Table 1: The OPAL 95% CL upper limit on the probability P for g frag- 
mentation as a function of m~. 



As always, it is important to keep in mind that if the decays to a neutral 
state of any kind with a lifetime shorter than ~ 10"'' sec, then P is effectively 
zero since the R^ will decay before traversing the tracker. Assuming a sufficiently 
long lifetime for the i?^, the limits of Table |I| can be interpreted in the context of 
the model for P described earlier. For P = 1, 1/2 and 1/4, one excludes m~ = 
1 — 20 GeV, 1.2 — 16.6 GeV and 1.9 — 13.6, respectively. We have already seen that 
the OPAL jets plus missing momentum analysis excludes m~ = 3 — 25 GeV for any 
P value not too close to 1; for P ~ 1, the upper limit dechnes to ~ 23 GeV. Thus, 
the limits from our analysis of the OPAL jets plus missing momentum channel are 
nicely complementary to the OPAL heavily-ionizing track limits; they confirm one 
another for a substantial range of m~. 



^^This is not quite the correct procedure in cases where both gluinos are accepted; the appro- 
priate multipHcation factor per gluino in that case is P — 
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7 The CDF signal for a penetrating charged gluino 
hadron 



The strength of this signal depends on the model used for gluino interactions and 
upon details of the detector. CDF's central muon system consists of two muon 
detection scintillators separated by iron. To be identified as a penetrating charged 
particle, a particle must (a) penetrate the iron, (b) be charged at the scintillator 
layer just before it enters the iron and (c) be charged at the exit detection layer. 
To be identified as a heavily-ionizing particle, the particle must also be charged as 
it exits from the primary interaction and its ionization must be clearly larger than 
minimal. 

Let us recall the picture we shall employ for the gluino as it traverses the de- 
tector. As in the OPAL analysis, the primary produced g is assumed to have some 
probability P to fragment (immediately) to a charged i?^-hadron. The ionization 
of the will be measured shortly after emerging from the interaction vertex. The 

then undergoes a certain number of hadronic interactions as it passes through 
the calorimeters before arriving at the inner muon detection layer preceding the 
iron. As described earlier, we imagine that at each hadronic interaction the light 
quark's and/or gluons are stripped from the i?-hadron (whether neutral or charged 
at the time) leaving the bare gluino which then has the same probability P to again 
become charged. Thus, the probability that the i?-hadron is charged just before 
entering the muon iron is again P. As it traverses the iron it will undergo sev- 
eral more hadronic interactions and so the probability that it exits as a charged 
i?-hadron is once again P. Altogether, we must reduce the cross section (after 
cuts to be discussed below) by P^. Once again, this assumes that all the possible 
charged i?-hadron states are effectively stable as they travel through the detector. 
If they decay rapidly to the or another neutral state, then this must be taken 
into account by an appropriate reduction of P. 

Whatever the value of P, we compute the event acceptance efficiency as follows 
HI. For a given m~, we generate events using ISAJET. We impose the triggering 



requirement that at least one of the g^s has 

Ir^l < 0.6 andpT > 15 GeV. (19) 

An efficiency of 0.8 is included for triggering on such a g. We next demand that 
at least one of the g^s satisfy the following heavily-ionizing, stable charged particle 
'reconstruction' requirements: 

\ri\ <l.O, IpI > 35 GeV , /3 > /3mm 
/37 < 0.85 for m~ > 100 GeV or p-f < 0.7 for m~ < 100 GeV . (20) 

We note that the jSj < 0.7 requirement we impose for m~ < 100 GeV is such that 
only events in which ionization is at least three times minimal (as compared to 
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Figure 29: The effective cross section aefr for one or more g to pass the heavily- 
ionizing penetrating particle cuts of Eqs. (|19D and (^), including the efficiencies 
quoted in the text. 



twice minimal if only /?7 < 0.85 is required) are accepted. This cut is stronger 
than that of the actual analysis |^^. We do this in the hope that the background 
will be even smaller than the conservative number used later. In the above, we 
use Puiin as given by the solid curve in Fig. |^. For P substantially smaller than 1, 
this is quite conservative given that ionization energy loss will be much less than 
that employed in the figure, which is for P = 1. Also, because we use /?min for 
P = 1 and because typical (3 values are substantially above Pmin, this analysis is 
quite insensitive to the choices of At and (AE) case. Finally, an efficiency of 0.5 is 
included for the reconstruction. Note that one g could provide the trigger but fail 
the reconstruction while the other g could pass the reconstruction cuts. In Fig. |2D 
we plot the effective cross section aefr as a function of m~ after including the above 
cuts and efficiencies, but before including P^. We note that no events pass the 
cuts for m~ < 50 GeV; the cuts would have to be weakened, which might result in 
the introduction of substantial background. 

In Ref. [E5|, it is stated that there are zero background events in L = 90 pb~^ 
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Figure 30: The 95%CL upper limit on the probabihty P for a ^ to fragment to a 
singly- charged hadron after production and collision is given as a function of 
and 10 background events. 



m- for Nb 



of data after the mass > 100 GeV cuts. The background level probably increases 
gradually as one lowers the m~ value considered down to 50 GeV. (Current cuts 
do not allow sensitivity below this.) However, even for the less stringent (3'y < 0.85 
cut the background level is estimated at < 12 events |^ for m~ = 50 GeV. To 



illustrate the situation, let us consider the cases of A''^ = and 10 background 
events. At 95% CL we require LP^a^s < 3 (A^b = 0) or < 7 (Nb = 10). The 
resulting 95% CL upper limits on P are plotted as a function of m~ in Fig. ^ 



We see that the limits on P are significant. In particular, for m~ ~ 50 GeV and 



10, we find that P > 0.09 is excluded. For m~ ~ 100 GeV and Nb = 0, 



P > 0.1 is ruled out, rising to P > 0.2 for m~ ~ 150 GeV. For m~ > 50 GeV, 
this result confirms the RunI jets-t-jr/r analyses that exclude values of m~ below 
130 - 150 GeV down to < 20 GeV for any P for SCI At and {AE) case choices. 
The heavily-ionizing track signal improves (though only slightly) for SC3 choices, 
and thus excludes m~ > 50 GeV (up to very big values) for P > 1/2 {i.e. for P 
values such that the jets+^T signal fails for the SC3 choices). We expect that, at 
large P, a CDF heavily-ionizing track analysis with weakened cuts would probably 
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be able to extend the excluded m~ range down to the OPAL m~ ~ 22 — 25 GeV 
lower bound (that applies for any P) based on the OPAL jets-f^y analysis and 
probably also down to the ~ 20 GeV bound (that applies for large P) from the 
OPAL heavily-ionizing track search. In any case, currently the only significant 
window for a ^-LSP in the P-m~ parameter space arises for SC3 choices and 
P > 3/4. The window at P ~ 3/4 is 25 GeV < m~ < 50 GeV, widening to 
23 GeV < m~ < 50 GeV for P ~ 1. 



8 A gluino NLSP decaying to gluon plus grav- 
itino 

For completeness, we consider the scenario in which the gluino is not the LSP, but 
rather the NLSP (next-to-lightest supersymmetric particle), with the gravitino 
(G) being the (now invisible) LSP. Such a situation can arise in GMSB models, 
including that of Ref. 0. In this scenario, the gluino decays via g gG. Early 
universe/rare isotope limits are then irrelevant. Further, the decay will be prompt 
from the detector point of view if is in the < few eV region such that the G 
is guaranteed to have no impact on Qh'^ (If the scale of supersymmetry 



breaking is so large that the g — gG decay lifetime is long enough that most ^'s exit 
the detector before decaying, then the results of previous sections apply.) The first 
examination of this scenario at the Tevatron appears in Ref. We are unaware 



of any studies of this scenario for the qqgg final state at LEP or LEP2. Here, we 
will give the 95% CL excluded mass domains based on the previously considered 



jets + missing momentum analyses of OPAL |3^ and CDF ||36|, ^7j. In our analysis, 
we will assume that the branching ratio of ^ — > gG is 100% (as appropriate if the 
g is the NLSP), and that the decay is prompt. We will also assume that the G 
has negligible mass compared to m~, and that other supersymmetric particles are 
much heavier than the gluino. 

Consider first, the OPAL analysis. Using exactly the same procedures and cuts 
as discussed earher in section 4, but applied to e~^e~ ^ Z — » qqgg qqgg +i/t, 
we have determined the efficiency for event acceptance and the resulting 95% CL 
upper limit on BR{Z qqgg)- These results appear in Fig. ^ Gluino masses 
below about 26 GeV are clearly excluded. 

For our CDF-based analysis of the g —>■ gG scenario we employ the same 
procedures as in section 5. We compute jets+j/r rates based on pp — > ggX. The 
plots analogous to those given earlier appear in Fig. |3^. We observe that the 
iets-^T signal cross section (after cuts) for a ^-NLSP is even larger than in the 
reference MSSM model. All values of m~ < 240 GeV (down to very small values 
that clearly overlap the OPAL exclusion region for this scenario) can be excluded 
at 95% CL based on the CDF L = 19 pb~^ data sample analysis. This result is 
stronger than the bound obtained in Ref. ESI . The same CDF analysis procedures 
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Figure 31: In the upper window, we plot the OPAL qTjgg event efficiency (after 
all cuts) in the g — > gG scenario. The lower window gives, as a function of m~, 
the corresponding 95% CL upper limits compared to the theoretical prediction for 
BR{Z^qqgg). 
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Figure 32: The cross section (after cuts) in the jets+]/r channel is compared to (a) 
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as a function of m~ for the g —>■ gG case. 
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applied at Runll will be able to exclude m~ values up to about 280 GcV. Analyses 
optimized for such higher masses will presumably be able to do even better. 

Overall, it is clear that a gluino NLSP decaying to gluon plus light gravitino 
can be excluded for essentially all m~ < 240 GeV. 

9 Insights for other new physics analyses 

In this section, we wish to emphasize a few interesting possibilities for other anal- 
yses for new physics that can be extracted from the lessons learned in our specific 
studies. 

The primary point to note is that our results imply that the jets plus missing 
momentum signal is immediately applicable for pair production of any type of 
stable or semi-stable {i.e. stable within the detector) neutral or charged heavy 
particle that is produced via the strong interactions. Examples of such particles 
abound in the literature. 

• Gauge-mediated supersymmetry breaking models can contain colored mes- 
sengers in the gauge-mediation sector that are stable or semi-stable. 

• In models with extra generations, one or more of the heavy quarks could be 
long-lived. 

• Semi-stable, strongly-interacting massive particles are proposed as a source 
of ultra-high-energy cosmic ray events. 

Pair production of a heavy stable particle produced via strong interactions gives 

rise to a substantial missing momentum signal due to the mismatch between the 
true momentum of each produced particle and the apparent energy of the jet 
associated with the particle (as measured after including calorimeter response and 
possible identification of any associated charged hadron track as a muon track 
within the jet). Further, the net missing momentum in a typical pair-production 
event does not tend to be aligned with the visible energy of the jet associated 
with any one of the heavy particles. This is because in a realistic Monte Carlo the 
pair-production process initiated by quarks and/or gluons in the colliding hadrons 
is accompanied by additional jets with high transverse momentum coming from 
initial state 'radiation'. 

In the case of pair production of heavy stable quarks at a hadron collider, the 
limits from the jets plus missing momentum analysis would be very complementary 
to the heavily-ionizing, penetrating track limits that rely more heavily on substan- 
tial modeling of the charge exchange and fragmentation for a heavy quark as it 
passes through the detector. As discussed earher, the rate for the latter signals 
scales roughly as P^, where P is the probability for the heavy quark to fragment to 
a charged (as opposed to neutral) heavy hadron. For small enough P, the missing 
momentum signal will be stronger than the penetrating track signal. In addition. 
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there is a very interesting hybrid signal that should be analyzed. A missing momen- 
tum trigger could be used to isolate events in which to look for a heavily-ionizing 
track.[^ This could be more efficient than the present CDF analysis which requires 
a penetrating track in order to have a trigger rate such that all events can be 
accepted. The jets plus missing momentum trigger would eliminate the need to 
require a penetrating track and one could just search for a heavily-ionizing track 
in events accepted by the trigger. The advantage would be that the probability for 
the heavily-ionizing track (without requiring penetration) scales only as P (rather 
than P3). 

It might be possible to take direct advantage of the mismatch between different 
ways of measuring the momentum of a heavy particle that is contained in a charged 
state after the initial interaction. The tracker would measure the true momentum 
of the particle. There are then two possibilities. 

• If the additional tracks are not present that cause the track observed in 
the tracker to be deemed as having penetrated to the muon detector, then 
this true momentum could be directly compared to the momentum of the 
particle as determined by the calorimeter response. We have seen that there is 
generally a very substantial difference. This situation would have probability 
oc P(l — P^) (including the probability for the initial track in the tracker). 

• Alternatively, if the track observed in the tracker is deemed to have pene- 
trated to the muon detector, one could compare the true momentum to that 
computed for the jet assuming the track belonged to a muon [see Eq. ([T5|) ]. 
The difference is substantial when the average (3 of the produced particle is 
large. 

In order to retain as many events as possible it would be best to use a simple multi- 
jet trigger (without necessarily requiring missing momentum). Of course, since we 
are looking for momentum discrepancies for a single jet, it would be necessary to 
perform a very careful study of backgrounds, such as that due to jets that are 
mismeasured and/or fragment to ii'^'s. 

10 Summary and Conclusions 

We have examined constraints on any model in which the gluino is the LSP. In 
section 2, we considered the relic cosmological density of a ^-LSP. We found that the 
relic density depends very strongly on the presence and nature of non-perturbative 
effects that could enter into the gluino and gluino-bound-state annihilation cross 
sections. Assuming a completely perturbative gg annihilation cross section leads 
to a relic density of Vth"^ ~ (m~/10 TeV)^. For m~ > 100 GeV, this level of 
relic density is probably inconsistent with bounds from limits from heavy isotopes, 

^^The g-LSP should also be searched for in the manner we describe. 
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underground detector interaction rates and the like. However, we found that non- 
perturbative effects can potentially decrease the relic density to flh^ ~ 10~^° for 
all m~ < 10 TeV, a level that would be entirely consistent with all constraints. 
Our conclusion is that, until the non-perturbative physics associated with gluino- 
gluino annihilation can be clarified, no reliable limits on the ^-LSP can be obtained 
from constraints requiring knowledge of its relic density. Thus, direct limits from 
accelerator experiments are of great interest. 

In section 3, we studied the manner in which a (stable) ^-LSP is manifested 
in a typical detector. The critical issue for experimental analyses is the average 
amount of visible momentum assigned to a gluino jet. For a given detector, this 
depends upon many ingredients, including the average hadronic collision length of 
the i?-hadron into which the g fragments, the average hadronic energy deposited 
in the various collisions experienced by the i?-hadron as it passes through the de- 
tector, and the typical velocity and charge of the i?-hadron. The hadronic collision 
length was estimated using the two-gluon model for total cross sections; one finds 
a collision length that is somewhat longer than for a typical light hadron. Colli- 
sion lengths that are twice as large and one-half as large as our central prediction 
were also considered. Two cross section models were employed for computing the 
average energy deposit (as a function of velocity) in each hadronic collision. The 
(generally fiuctuating) charge of the i?-hadron as it passes through the detector is 
also a crucial ingredient and is characterized in terms of the probability P for the 
g to turn into a stable charged R"^, such as gud, as opposed to a neutral state, 
such as the = gg, after a hadronic collision. Simple quark counting models 
suggest P < 1/2 and probably much smaller if the gg bound state is important. 
For P — 0, the energy (^momentum) assigned to a gluino jet will be equal to 
the amount of the g kinetic energy that is deposited in the calorimeters due to 
hadronic collisions. For P > 0, the ionization energy deposits must be included 
and the possible interpretation of an R^ track in the central tracker as a muon 
within the ^-jet must be taken into account. 

In order to do this properly in a Monte Carlo context, for any given value of 
P, the momentum measured for each g is computed on an event-by-event basis, 
including (for P 7^ 0, 1) random changes (according to the value of P) of the 
charge of the i?-hadron at each hadronic collision as the g passes through the 
detector. Procedures are highly dependent upon the detector and specific analysis 
in question. For example, in the LEP OPAL jets -|- missing momentum analysis, 
if the i?-hadron is an R^ in the tracker and penetrates as an R^ to the muon 
chamber, then the ^-jet is declared to contain a muon and a procedure for adding 
in the supposed muon track momentum (and correcting for its presumed minimal 
ionization energy deposit in the calorimeter) is followed. In contrast, in the CDF 
jets -|- missing momentum analysis for Tevatron RunI, if the i?-hadron is an 
in the tracker and appears as an in one of the muon chambers, and if the 
net measured calorimeter energy is not too large, then the ^-jet is declared to be 
muonic and the event is discarded. 
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We studied the momentum typically assigned to the g-^et as a function of P, 
for the g masses and velocities of relevance, in the OPAL and CDF analyses. For 
all P (for P < 1/2), we found that the CDF (OPAL) procedure implies that the 
momentum assigned to the ^-jet is (on average) only a small fraction of its actual 
momentum unless m~ is smaller than a few GcV. This is true even for the cross 
section choice that overestimates energy deposits and even though, in the OPAL 
procedure, we allow for the appropriate fraction of cases (determined by P) in 
which the g penetrates to the muon chamber and has an track that is treated 
as a muon component of the jet in reconstructing the jet energy. Thus, when the g 
is the lightest supersymmetric particle, the jets plus missing momentum signature 
at colliders is, indeed, relevant. In fact, this would be the dominant standard 
SUSY signal if all other supersymmetric particles, in particular those with strong 
production cross sections, are significantly heavier than the g. 

Section 3 ended with a discussion of the effects of incomplete containment 
of a shower from a hadronic interaction that takes place near the outer edge of 
the hadronic calorimeter (or outer edge of uninstrumented iron). Effects, on the 
OPAL and CDF analyses summarized below, from the failure to include the shower 
energy in the measured jet energy and from the extra tracks in the subsequent 
muon-chamber(s) are outlined. 

As noted, existing jets plus missing momentum analyses at both LEP and the 
Tevatron are relevant to excluding a range of m~ values in the ^f-LSP scenario. 
In section 4, we demonstrated that the OPAL LEP data analysis that has been 
performed in order to search for Z — > X1X2 (with X2 ~^ QQXi) i^i the jets plus missing 
momentum channel can be applied to Z — > qqgg events. For P = 0, 1/4, 1/2, 3/4, 
we found that m~ values from ~ 3 GeV up to ~ 25 GeV are excluded at the 95% CL, 
for all choices of path length Xt and (AP) energy loss (per hadronic collision) case 
considered. For P = 1, and after including energy smearing and fragmentation 
effects, the upper limit of the excluded range declines to m~ ~ 23 GeV for our 
standard or "SCI" choices of and (AP) case. There is almost no change of the 
excluded range of m~ for possible extreme choices of and (AP) (with scenario 
labels "SC2" and "SC3"). For the "SCI" choices, results for P ~ 1 are sensitive 
to whether or not we include energy smearing and fragmentation effects. If these 
effects are not included, the fluctuations in measured jet energy are reduced and 
no limit is possible for P = 1 from OPAL jets + missing momentum data. (But, 
as discussed below, much the same range of 3 < m~ < 20 GeV is excluded by the 
heavily-ionizing track signal.) In contrast, for P < 3/4, the excluded range of m~ 
is essentially independent of whether or not energy smearing and fragmentation 
are included. Turning to LEP2, we noted that accumulated luminosities will not 
be adequate to improve the LEP Z-po\e limits. A next linear collider operating at 
^/s — 500 GeV would be able to extend the LEP limits, but probably not beyond 
the limits that are imposed by our Tevatron analysis. 

In section 5, we analyzed constraints from the Tevatron, assuming that all 
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other SUSY particles are much heavier. We beheve the resulting limits on m~ to 
be conservative. We examined the jets plus missing momentum channel using cuts 
and procedures based on the currently pubhshed CDF analysis of L = 19 pb~^ 
of RunI data. The cross section hmits obtained by CDF translate to a range of 
excluded m~ values. At 95% CL, wc exclude m~ up to ~ 130 — 150 GeV (the precise 
upper limit depending on P) down to at least 20 GeV (at a very high CL), for 
"SCI" or "SC2" choices of Ay and {AE) case. For "SC3" choices (corresponding 
to long path length and small hadronic energy deposits per collision for the g) the 
current CDF analysis can only exclude the above range of m~ for P < 1/2. Thus, 
for all but "SC3" choices, the CDF RunI limit overlaps the OPAL limit for any 
value of P, and all values of m~ in the ~ 3 GcV— 130 GcV range are excluded. For 
"SC3" Xt and (AE) case choices, these same CDF limits apply only for P < 1/2. 
This lack of sensitivity of the CDF analysis at large P to long At and/or small 
(AE) could be ehminated by a re-analysis of the data that retains muonic jets. 

Runll Tevatron data in the jets plus missing momentum channel can be ex- 
pected to extend the exclusion region to higher masses; depending upon P, we 
found that roughly m~ < 160 - 180 GeV will be excluded for "SCI" or "SC2" 
choices of At and (AE) case. For "SC3" choices and high P, only m~ < 130 GeV 
would be excluded. Such sensitivity is substantially worse than that found for 
the MSSM with mSUGRA boundary conditions, for which one can probe out to 
roughly m~ < 250 GcV. Possibly the Runll reach in the ^'-LSP scenario could 
be extended if systematic errors are smaller than anticipated. The above limita- 
tion assumes that S/B > 0.2 is required for a detectable signal. Alternative cuts, 
with smaller B at high m~, might also yield a larger reach. Although we have not 
specifically performed the analysis, the Tevatron results suggest that the LHC can 
be expected to rule out a ^-LSP with m~ up to at least 1 TeV. 

We also explored limits on a ^-LSP deriving from the non-observation of a 
pseudo-stable charged track which is heavily-ionizing. The strength of such signals 
depends on P. In section 6, we reviewed the OPAL results. OPAL performed a 
direct search for such states using cuts in the dE/dx — |p| plane, concluding that 
for P ~ 1/2 (P ~ 1) one can exclude m~ in the ~ 1 — 17 GcV (~ 1 — 20 GeV) mass 
range. For heavy-ionization signals at higher masses we must turn to the Tevatron. 
CDF looks for events containing a pseudo-stable penetrating charged track which 
is heavily-ionizing. In section 7, we demonstrated that, depending upon P, g- 
pair production can lead to a significant cross section (after imposing the CDF 
cut, penetration and ionization requirements for identifying such events with small 
background). We have estimated the upper limit from RunI data on the probability 
P of charged fragmentation of the g. The upper limit can be roughly parameterized 
as P ~ 0.3(m~/200 GeV) for 100 <, m~ <, 250 GeV. For m~ < 140 GeV, this 
means that P < 0.18 is required. Meanwhile: the jets plus missing momentum 
limits based on OPAL and CDF analyses exclude 3 GeV < m~ < 130 — 150 GeV 
for P < 1/2; the OPAL jets plus missing momentum analysis excludes ~ 3 GeV < 
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m~ < 25 GeV for any P not too near 1 (~ 3 GeV < m~ < 23 GeV for P = 1); 
and, the CDF jets plus missing momentum analysis excludes m~ from ~ 20 GeV 
to ~ 130 GeV for P = 3/4 and P = 1 for all but "SC3" choices of At and {AE) 
case. For P > 1/2 (independent of At and (AE)), the CDF heavily- ionizing track 
analysis excludes 50 < m~ < 200 GeV. This leaves only the possibility that "SC3' 
choices apply, that P lies in the (less likely) P > 3/4 range, and that m~ lies in the 
~ 23 — 50 GeV window. Very probably, an extension of the CDF heavily-ionizing 
penetrating particle analysis with weakened cuts appropriate to these lower masses 
could exclude this window. 

For completeness, in section 8 we also considered the scenario where the gluino is 
the NLSP (next-to-lightest supersymmetric particle) and the gravitino is the LSP. 
Such a situation is quite possible in models with gauge-mediated supersymmetry 
breaking. In this scenario, the gluino decays via g gG and the G is invisible. 
There is then a strong jets-fj/y signal at both LEP and the Tevatron. We repeated 
the LEP OPAL-based analysis and the L = 19 pb^^ CDF-based analysis for this 
case and found that m~ < 240 GeV can be excluded at 95% CL. Runll should be 
able to extend the excluded region to at least m~ ~ 280 GeV. 

Finally, we urge our experimental colleagues to take note of our remarks in 
section 9 regarding the applicability of our procedures in the jets plus missing 
momentum channel, or hybrid procedures such as combining a jets plus missing 
momentum trigger with a heavily-ionizing track requirement, to placing limits on 
other exotic particles, such as a heavy stable quark. We also note that a search for 
heavily-ionizing tracks in events with jets plus missing momentum should prove 
very valuable for excluding P > 1/2 ^-LSP scenarios. 
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